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Abstract. In this lecture a short introduction is given into the theory of 
the Feynman path integral in quantum mechanics. The general formulation 
in Riemann spaces will be given based on the Weyl- ordering prescription, 
respectively product ordering prescription, in the quantum Hamiltonian. Also, 
the theory of space-time transformations and separation of variables will be 
outlined. As elementary examples I discuss the usual harmonic oscillator, the 
radial harmonic oscillator, and the Coulomb potential. 
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I Introduction 

It was Feynman's (and Dirac's [20]) genius [32, 33] to realize that the integral kernel 
(propgator) of the time-evolution operator can be expressed as a sum over all possible 
paths connecting the points q' and q" with weight factor exp [iS(q", q'; T)/h] , where 
S is the action, i.e. 

K(q",q';T)= ^ A jS{ q ", q ',T)/h (L1) 

all paths 

with some appropriate normalization A. 

Surprisingly enough, the same calculus ("same" in the sense of a naive analytical 
continuation) was already know to mathematicians due to Wiener in the study of 
stochastic processes. This calculus in functional space ("Wiener measure") attracted 
several mathematicians, including Kac (who mentioned being influenced by Feynman's 
work!), and was further developed by several authors, where best known is the work 
of Cameron and Martin. The standard reference concerning these achievements is the 
review paper of Gelfand and Yaglom [37] , where all these early work was first critically 
discussed. 

Unfortunately, the discussion between physicists and mathematicians remains 
near to nothing for quite a long time, except for [37]; the situation changed with 
Nelson [80], and nowadays there are many attempts to understand the path integral 
mathematically despite its pathalogies of "infinite measure" , "infinite sums of phases" 
with unit absolute values etc. 

In particular, the work of Morette-DeWitt starting with her early paper [76] gave 
rise to a beautiful theory of the semiclassical expansion in powers of h [18, 19, 77-79]. 
As it is known for quite a long time, the propagator can be expressed semi-classically as 
e iS C i/h^ w ith s cl the classical action, times a prefactor. This prefactor is remarkably 
simple, namely one has 



K WKB ( X ",x';t",t>) = [g( X ')g( X ")]-^-l 



D/2 

h 



g = det (gab) of some possible metric structure of a Riemannian space, and the 
determinant 

, ( d 2 Sci[x',x"}\ . . 

M ■= det ( - a ? L ' (1-3) 



v dx' a dx' b j 

is known as the Pauli-van Vleck-Morette determinant. The semiclassical (WKB-) solu- 
tion of the Feynman kernel (we use the notions semiclassical and WKB simultaneously) 
is based on the fact that the harmonic oscillator, respectively, the general quadratic 
Lagrangian, is exactly solvable and its solution is only determined by the classical 
path and not on the summation over all paths. As it turns out, an arbitrary kernel 
can be expanded in terms of the classical paths as an expansion in powers of h. The 
semiclassical kernel (at least its short time representation) is known since van Vleck, 
derived by the correspondence principle. Later on Pauli [83] has given a detailed dis- 
cussion in his well-known lecture notes. More rigorously the short time propagator for 
the one-dimensional case was discussed by Morette [76] in 1951 and a few years later 
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by DeWitt [17] extending the previous work to D-dimensional curved spaces. Con- 
cerning the semi-classical expansion, every path integral with a Hamiltonian which 
is quadratic in its momenta can be expanded about the semiclassical approximation 
(1.2) giving a consistent and converging theory, compare DeWitt-Morette [18, 19, 74, 
75, 77]. 

Supersymmetric quantum mechanics provides a very convenient way of classifying 
exactly solvable models in usual quantum mechanics, and a systematic way of 
addressing the problem of finding all exactly solvable potentials [56] . 

By e.g. Dutt et al. [16, 29] it was shown that there are a total of twelve different po- 
tentials. A glance on these potentials shows that their corresponding Schrodinger equa- 
tion leads either to the differential equation of the confluent hypergeometric differential 
equation with eigen-functions proportional to Laguerre polynomials (bound states) 
and Whittaker functions (continuous states), respectively to the differential equation 
of the hypergeometric differential equation with eigen-functions proportional to Ja- 
cobi polynomials (bound states) and hypergeometric functions (continuous states). 
In Reference [16] this topic was nicely addressed and it was shown that in principle 
the radial harmonic oscillator and the (modified) Poschl- Teller-potential solutions are 
sufficient to give the solutions of all remaining ones, together with the technique of 
space-time transformations as introduced by Duru and Kleinert [27, 28, 62] (see also 
reference [49], Ho and Inomata [55], Steiner [90], and Pak and Sokmen [82]) enables 
one to give the explicit path integral solution of the Coulomb V^ c '{r) = —e 2 /r (r > 0) 
and the Morse potential V^ M \x) = (H 2 A 2 /2m)(e 2x — 2ae x ) (x G R) (compare e.g. 
reference [46] for a review of some recent results) . The same line of reasoning is true for 
the path integral solution of the (modified) Poschl- Teller potential [5-35] which give in 
turn the path integral solutions for the Rosen-Morse V^ RM ^ (x) = A tanhx— Bj cosh 2 x 
(i 6 R), the hyperbolic Manning- Rosen V^ MRa \r) = — Acothr + Bj sinh 2 r (r > 0), 
and the trigonometric Manning- Rosen-like potential V^ MRb ^ (x) = —A cot x + B/ sin 2 x 
(0 < x < 7r), respectively [44]. 

These lecture notes are far from being a comprehensive introduction into the 
whole topic of path integrals, in particular if field theory is concerned. As old as 
they be, the books of Feynman and Hibbs [34] and Schulman [86] as still a must for 
becoming familiar with the subject. A more recent contribution is due to Kleinert 
[64]. Myself and F. Steiner are presently preparing extended lecture notes "Feynman 
Path Integrals" and a "Table of Feynman Path Integrals" [50, 51], which will appear 
next year. 

Several reviews have been written about path integrals, let me note Gelfand and 
Jaglom [37], Albeverio et al. [1-3], DeWitt-Morette et al. [19, 79], Marinov [73], and 
e.g. for the topic of path integrals for Coulomb potentials [46] . 

The contents of the lecture is as follows. In the next Chapter I outline the basic 
theory of the Feynman path integral, i.e. the lattice definition according to the Weyl- 
ordering prescription in the Hamiltonian and a related prescription which is of use 
in several applications in my own work. Furthermore, the technique of canonical 
coordinate-, time- and space-time transformations will be presented. The Chapter 
closes with the discussion of separation of variables in path integrals. 

In the third chapter I present some important examples of exact path integral 
evaluations. This includes, of course, the harmonic oscillator in its simplest form. I also 
discuss the radial path integral with its application of the exact treatment of the radial 
(time-dependent) harmonic oscillator. The chapter concludes with a comprehensive 
discussion of the Coulomb potential, including the genuine Coulomb problem, its D- 
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dimensional generalization, and a generalization with additional axially-symmetric 
terms. 
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1. The Feynman Path Integral 

In order to set up the requirements of the path integral formalism we start with the 
generic case, where the time dependent Schrodinger equation in some D-dimensional 
Riemannian manifold M with metric g ab and line element ds 2 = g ab dq a dq b is given by 

n 2 . TW J T , s h d 



%i) = TS *tef). (11) 



\1/ is some state function, defined in the Hilbert space C 2 - the space of all square 
integrable functions in the sense of the scalar product (/i, f-2) = J M y/g fi(q)f2(q)dq 
[g := det(^ab), /i, f 2 G £ 2 ] and A^s is the Laplace-Beltrami operator 



A LB := g-td a gtg ab d b = g ab d a d b + g ab (d a In ^)d b + g ab a d b (1.2) 

(implicit sums over repeated indices are understood). 

The Hamiltonian H := — ^-Alb + V(q) is usually defined in some dense subset 
D(H) C £ 2 , so that if is selfadjoint. In contrast to the time independent Schrodinger 
equation, H^> = E~^>, which is an eigenvalue problem, and equation (1.1) which are 
both defined on D(H), the unitary operator U{T) := e - lTH / h describes the time 
evolution of arbitrary states \1/ G C 2 (time-evolution operator): H is the infinitessimal 
generator of U. The time evolution for some state \1/ reads: ^{t") = U(t",t')^f(t') 
Rewriting the time evolution with U (T) as an integral operator we get 

*(«"; t") = J VgW)K(q", q'; t", t'Mq'; t')dq> \ (1.3) 

where K(T) is the celebrated Feynman kernel. Equations (1.1) and (1.3) are 
connected. Having an explicit expression for K(T) in (1.3) one can derive in the 
limit T = e — > equation (1.1). This, on the other hand proves that K{T) is indeed 
the correct integral kernel corresponding to U (T). A rigorous proof includes, of course, 
the check of the selfadjointness of H, i.e. H = H* . Due to the semigroup property 

U(t 2 )U(t l )=U(t l +t 2 ) (1.4) 

(H time-independent) we have 

K(g", q'; t" + t') = J K(q", q; t" + t)K(q, q'; t + t')dq (1.5) 

In the case of an Euclidean space, where g ab = 5 ab , S is just the classical action, 
S cl = J [f q 2-V(q)]dt = J C ci {q,q)dt, and we get explicitly [AqW := (q^-q^-V), 

qti) = q (tj), tj = t' +je, e = (*" - t')/N, N -> oo, D = dimension of the Euclidean 
space] : 



ND AT 1 




^A 2 qW -eV(qW) 
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«(*")=«" 



/ ^HU 



m 



q 2 - V(q) 



dt 



(1.6) 



q(t')=q' 



For an arbitrary metric g ab things are unfortunately not so easy. The first formulation 
for this case is due to DeWitt [17]. His result reads: 

K(q",q';T) 

9(t")=q" , _ 

= J y/gV DeW q(t)exp Ji J 

q{t>)=q' ^ 1 

/ m \ ND/2 ^il} r I 



rn 
~2 



(j) 



x 



i N 

exp u£ 



^g ab (q(i-V)Aq a ^Aq b ^ - eV(q^) + e^R(q^) 
2e bm 



(1.7) 



(i* = <7 ab (r^, c - r^ >a + 1^ - T d cb T c ad ) - scalar curvature; Y a hc = g ad (g bd , c + g dc , b - 

gbc,d) - Christoffel symbols). Of course, we identify q' = q^ and q" = q( N ) in the 
limit N — > oo. Two comments are in order: 

1) Equation (1.7) has the form of equation (1.6) but the corresponding S = j Cdt 
is not the classical action, respectively the Lagrangian C is not the classical 



Lagrangian C(q,q) = Y9abQ a Q ~ ^(<z)j but rather an effective one: 



'e// 



= J C e ffdt = J (C c i - AV DeW )dt. 



(1.8) 



The quantum correction AVoeW = — gm-R is indispensible in order to derive from 
the time evolution equation (1.3) the Schrodinger equation (1.1) (e.g. [24]). The 
appearance of a quantum correction AV is a very general feature for path integrals 
defined on curved manifolds; but, of course, AV ~ h 2 depends on the lattice 
definition. 

2) A specific lattice definition has been chosen. The metric terms in the action are 
evaluated at the "prepoint" q^- 1 ) . Changing the lattice definition, i.e. evaluation 
of the metric terms at other points, e.g. the "postpoint" q^ or the "mid-point" 
q(i' := \{q^' + q^~ 1 ^) changes AV, because in a Taylor expansion of the relevant 
terms, all terms of 0(e) contribute to the path integral. This fact is particularly 
important in the expansion of the kinetic term in the Lagrangian, where we have 
A*qV)/ e ~0(e). 

One of the basics requirements of the construction of the path integral is Trotter's 
[93] product formula. Discussions and proofs can be found in many textbooks on 
functional integration and functional analysis (e.g. Reed and Simon [85], Simon [87]). 
Let us shortly notice a simple proof [87] : 

Theorem: Let A and B be selfadjoint operators on a separable Hilbert space so that 
A + B, defined on D(A) . . . D(B), is selfadjoint. Then 



e iT(A+B) = s _ Um ( e itA/N e itB/N\ 



N 



(1.9) 
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If furthermore, A and B are bounded from below, then 



-t(A+B) 



N 



= s-\im(e- tA / N e- tB / N ) , 



(1.10) 



Proof ([80,87]): Let S T = e iT(A+B \ V T = e iTA , W T = e iTS , U T = V T W T , and 
let \&t = St^ for some $ G with underlying Hilbert space Ti. Then 



\\{ST-u» /N )n = 



N-l 



J2u^ /N (s T/N -u T/N ) N -^ 



3=0 



< N sup \\(S T - U T/N )y\\. 

0<s<T 



(1.11) 



Let $ G D(A) n Then s" 1 ^ - 1)$ -> K A + B )® for s and 



-> iS$ + iA$ + 



(1.12) 



hence 



lim 



iV||(5 T/ jv - U T/N )Q\\ -> 0, for each $ e D(A) n (1.13) 



Now let £) denote D(A) n with the norm ||(A + + ||$|| = ||$|U+s- 

By hypothesis, D is a Banach space. With the calculations of equations (1.11- 
1.13), {N(S t / N — U t / N )} is a family of bounded operators from D to H with 
sup JV {A^||(5' T / A r — U T / N )&\\} < oo for each $. As a result, the uniform boundedness 
principle implies that 



N\\(S T /N-U T/N )$\\<C\m A + B 



(1.14) 



for some positive C Therefore the limit (1.14) is uniform over compact subsets of 
D. Now let f G D. Then s — * \l/ s is a continuous map from [0, T] into .D, so that 
{^ s |0<s<T}is compact in D. Thus the right-hand side of equation (1.14) goes to 
zero as N — > oo The proof of equation (1.10) is similar. 
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The conditions on the operators A and B can be weakened with the requirement 
that they are fulfill some "positiveness" and possesses a particular self-adjoint exten- 
sion. 

Thus we see that ordering prescriptions in the quantum Hamiltonian and corre- 
sponding lattice formulations in the path integral are closely related (compare also 
[22, 23]). This can be formulated in a systematic way (e.g. [68] ): 
We consider a monomial in coordinates and momenta (classical quantities) 



M(n,m) =q» 1 ...q^p Ul . . .p v " 
We want to have a correspondence rule according to 



exp 



I v 



Dq(u, v; q, p) = Q(u, v) exp 



-(uq + vp) 



(1.15) 



(1.16) 



to generate operators q, p from coordinates q, p. This produces the mapping 



M(n,n) 



d n+m D n (u,v;q,p) 



<9w Ml . ..du^dv 



(1.17) 



u=v=0 



Some better known examples are displayed in the following table 



Correspondence Rule 



0(w, v) 



Ordering Rule 



Weyl 
Symmetric 
Standard 
Anti-Standard 
Born- Jordan 



u ■ v 



cos 



1=0 v 7 

i(qV-+p"V) 



exp — i 



u ■ V 



exp I i 



2 

u ■ V 



q n p m 



p m q n 



U ■ V u ■ V 



sin 



2/2 



m + 



1=0 



We can make this correspondence explicit. Let us consider the Fourier integral 



A(p, q) = J A(u, v) exp 



I v 



dudv 



A(u,v) 



We define the operator A n (p, q) via 



A(p, q) exp 



- — {uq + i v P) 



(1.18) 



dpdq 



-4 n (p,q)= J M u > v)Q(u, v) exp 



-(wq + ivp) 



dudv, 



(1.19) 



II.2 



Weyl-Ordering 



which gives 
A"(p,q) = 



A(p, q)Q(u, v) exp 



-^u(q-q)-^v(p-p) 



(2tv%) d 

The inverse transformation is denoted by A & {p, q) and has the form 



dudvdpdq. (1.20) 



^u(q-q) + i^v(p-p) 



In particular, the 5-function has the correspondence operator 

1 



A n (p-p,Q-q) 



fl(u, v) exp 



- ^u(q-q) - ^v(p-p) 



dudv. (1.21) 



dudv. (1.22) 



(2irh) D 

Let us turn to the Hamiltonian. Given a classical Hamiltonian H c i (p, q) the 
quantum Hamiltonian is calculated as 



H(p, q) 
Hd(u,v) = 



exp ( —up + —vq ) Q(u, v)H c i(u, v)dudv 



1 



exp 



^up - ^uq )H c i(p, q)dpdq. 



(1.23) 



(2nh) 2D 

Obviously, by proposing a particular classical Hamiltonian depending on variables q' 
and q", respectively, produces a particular Hamiltonian function, 

H(p,q'\q') = j^ E ^ J il{q" -q\y)^y/ h H cl [p,W + q')-v]dudv (1.24) 

For Q = 1 and O = cos respectively, we obtain Hamiltonian functions according to 

H(p,q",q')=H(p, 1 1 (q' + q")) 
H(p,q",q') = ±[H c i(p,q") + H c i(p,q')) 
which are the matrix elements 

<</|B(p,q)|</'> (1.26) 

of the quantum Hamiltonians which are Weyl- and symmetrically ordered, respectively. 
Choosing Cl(u,v) = exp[i(l — 2u)uv] [54] yields 

H(p, q; u) = £-g ab (q)p aPb + - u)p a g a \ h ~^{\- «) V% + AV Weyl . (1.27) 



(1.25) 



2m ' m 
with the well-defined quantum potential 



Weyl 



(9 ab r d ac n d - R) 



g ab T a T b + 2(g ab T a ), b + g a » ab 



ab-i 



(1.28) 



8m w ac ua ' 8m 
u= \ corresponds to the Weyl prescription and is clearly emphazised. 

In the next two Sections we discuss two particular ordering prescriptions and the 
corresponding matrix elements which will appear as the relevant quantities to be used 
in path integrals on curved spaces, it will be the Weyl-ordering prescription (leading 
to the midpoint rule in the path integral) and a product ordering (leading to a product 
rule in the path integral). 

Another approach due to Kleinert [63, 64] I will not discuss here. 



9 



General Theory 



2. Weyl-Ordering 

A very convenient lattice prescription is the mid-point definition, which is connected to 
the Weyl-ordering prescription in the Hamiltonian H. Let us discuss this prescription 
in some detail. First we have to construct momentum operators [83]: 

which are hermitean with respect to the scalar product (Ji,/^) = / f\ft\f9^1- ^ n 
terms of the momentum operators (2.1) we rewrite the quantum Hamiltonian H by 
using the Weyl-ordering prescription [49,69,72] 

H(p, q) = ^-(g ab PaPb + 2p a g ab Pb+PaPbg ab ) + AV Weyl (q) + V(q). (2.2) 
Here a well-defined quantum correction appears which is given by [49,72,81]: 

Wweyi = ^(9 ab r d ac r c bd -R) = ^ [9 ab r a r b + 2(^r a ), b + g ab ab ] (2.3) 

The Weyl-ordering prescription is the most discussed ordering prescription in the 
literature. Let us start by defining it for powers of position- and momentum- operators 
q and p, respectively [69]: 

(1 \ m _ m l 
5) E^^9"-W. < 2 - 4 > 

with the matrix elements (one-dimensional case) 

(1 \ m i 
2) E , !(m "l iy 1 m -' <</'\p'W>l' 

and all coordinate dependent quantities turn out to be evaluated at mid-points. Here 
was used that the matrix elements of the position \q > and momentum eigen-states 
\p > have the property 

< q"W >= (g'g'T 1/4 $(q" - ?'). < q\p >= {^h)~ D/2 e ipq/h . (2.6) 

This power rule is nothing but a special case of a more general prescription. Of course, 

well behaved operator valued functions are included. 

We have the Weyl-transform of an operator A(p, q) as [72]: 

A(p,q) = J dvjv v <q-l\A{v,q)\q+l >, (2.7) 

10 



II. 2 Weyl-Ordering 
with the inverse transformation 

^ = (2nh) D J dpdqA & ti A (p> 9) 

A(p,q) = J due ipu |p — f ><p+f|. 
Let us consider some simple examples. First of all for some operator /(p) 

y cfc e [pv / h < q- f|/(p)k+ f > 

cfcc^' e ip " /a < g - f |/(q)b' >< p'|g + f > 

J dvdp'e^- p >/ h f(p') = f{p). (2.9) 

Hence (and similarly) 

/(p) <=► /(p), /(g) <=► /(«)• (2.10) 

Straightforwardly one shows the following correspondence 



/ 



1 fl 



(2.11) 

Special cases are 

P,F(q)P 3 ^ ^ + __j^ + __jjr (9 ). 
For the case that -F(g) is some symmetric N x N matrix we have 



p.p .i^(q) + 2p ; /•"' (<i)/^ + i^(q)p.p. 



PiPjF^iq). (2.13) 



There is, of course, a one-to-one correspondence between the function A(p, q) and the 
operator A(p, q), called Weyl-correspondence. It gives a prescription how Weyl- 
ordered operators can be constructed by the classical counterpart. Starting now 
with an operator A, the Weyl-correspondence gives an unique prescription for the 
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construction of the path integral. We have for the Feynman kernel for an arbitrary 
iV G N [which is due to the semi-group property of U(T), i.e. U(t 1 +t 2 ) = U(ti)U(t 2 )]: 



K(q",q';T) = (q 



exp 



h 



H(p,q) 



exp 



^H(p, 9 ) 



(2.14) 



Observing 



< g'|A(p,q)|g" > = J du e iuq < q'\p-% ><p + f \q" > 

du exp 



= e^'-^/ h 5(q- q,+ 2 q '' ) 



1 / / r,\ , 1 ( Q +Q 



(2.15) 



and making use of the Trotter formula e U ( A + B ) = s — liniTv^oo (e ltA / N e ltB / N ^ 
and the short-time approximation for the matrix element < q" | e 



9 >: 



< q^' | exp 
1 



ieH(p,q)/n 



> 



(27r^)° 
1 

(2yrn) 



,0) 



< 9 
, , , dp exp 



dpdq e lh (P^/ h J dudv exp [i(q — q)u + i(p — p)v] 



,0-i) 



> 



jpiq^-q^-jHip,^) 



le 



(2.16) 



where q^ = \{q^ + q^ 1 - ) ) is the mid-point coordinate. 
Let us consider the classical Hamilton-function 



Wci(p, 9) = ^9 ab (q) [Pa ~ A a (q)] [p b - AM] + V(q). (2.17) 

The corresponding quantum mechanical operator is not clearly defined due to the 
factor ordering ambiguity. Let us try the manifestly hermitean operator 

H(p, q) = ^-g-~* [p a - A a (q)]^(q)^ b (q) [p - A (q)]<r*(q) + V(q). (2.18) 



The Weyl-transformed of H reads 

H(p, q) = H c i(p, q) + ^ 



- g *{g*g ab ) 



,ab 



4(f i ).j>- i ), l »v'-!(s- i )>" 
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= n cl (p,q) + —[rr a r l mb g ab -R] 



= H c i(p, q) + £- g ab r a r b + 2(g ab r a ), b + g a \ ah 



8m 



abi 



-,ab 



(2.19) 



and, of course, AVweyi appears again. For the path integral all quantities have to 
be evaluated at q^K Note: The Weyl-transformed of the the Weyl-ordered operator 
(2.16) is just the classical Hamiltonian (2.17). Inserting all quantities in equation 
(2.13) we obtain the Hamiltonian path integral: 



N-l „ N , (j) 

K(q",q>;T) = WWT* ^ \{ / • \[ / * 

7 = 1 7 = 1 



' N 

x exp 



(2.20) 



The effective Hamiltonian to be used in the path integral (2.20) reads, 



H(p, q) = n eff (p^,q^) = l-g-\q^)p^p^ + V(q^) + AV Weyl (q^). (2.21) 
With the help of the famous Gaussian integral 

£j^**dx=^ew(J^) (2-22) 

and, respectively, by its D-dimensional generalization 

dp e^-y a "^ = (2yr) D / 2 A /det(^) exp ( - ^/«6«V) , (2.23) 

we get by integrating out the momenta the Lagrangian path integral which reads 
(MP = Mid-Point): 



q(t")=q" 



K(q",q>;T) = [g(q>)g(q")}-* J y/gV M pq(t)exp C eff (q,q)dt 



q(t')=q> 

ND I N-l 



: = ww<r>r* ^ 2 (II / 



X 



N , . r 

n yjgi^i)) exp \ i ^-g ab {q^)Aq a ^Aq b ^ - eV(q^) - eAV Wey i(q {j) ) 



2e 



(2.24) 
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Equation (2.24) is, of course, equivalent with equation (1.7). The mid-point pre- 
scription arises here in a very natural way, as a consequence of the Weyl-ordering 
prescription. This is a general feature that ordering prescriptions lead to specific 
lattices and that different lattices define different AV. 

To proof that the path integral (2.24) is indeed the correct one, one has to show 
that with the corresponding short time kernel 



2tt i eh 



x exp 



h 



m 
2e [ 



, ab (q^)Aq a ^Aq b ^ - eV(q^) - eAV Weyl (q^> 



(2.25) 



and the time evolution equation (1.3) the Schrodinger equation (1.1) follows. For this 
purpose a Taylor expansion has to performed in equation (1.3) yielding 



di 



B *(q";t) + B qb 



d^(q";t) 



q dq' h 



+ B„a 



d 2 V(q";t) 



11 X„ib,a 



dq'°q' c 



+ ..., (2.26) 



where the coefficients in the expansion are given by 



Bo = J dq'^)K{q",q'-e) 



2n i eh 



x / dq'g*(q)g*(q')exp ( ^-Aq a a g b Aq b 



i m 



B qb = I dq'yfgJq T )K(q",q>;e)Aq b 



{^£) D/2 9~H<l") I dq'gHq)gHq')e^(^Aq a a g b Aq b )Aq b 



(2.27) 



(2.28) 



B qaqb = I dq'yfgJ4)K{q",q'-e)Aq a Aq l 



<T 4 («") J dq'g^ (q)g* (q') exp Aq a a g b Ag b j Aq a Aq b ^ 



From these representations it is clear that we need equivalence relations (in the sense 
of path integrals) for Aq a Aq b etc. They are given by 



Aq a Aq b =—g ab 
m 



(2.30) 



Aq a Aq b Aq c Aq d =(^-^j [g ab g cd + g ac g bd + g ad g bc ] . (2.31) 
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Aq a Aq b Aq c Aq d Aq e Aq f 



i eh V 
m J 



9 ab 9 cd 9 cd + 9 ac 9 M 9 ef + 9 ad 9 bc 9 ef 



+ 9 ab 9 ce 9 df + 9 ab 9 cf 9 de + 9 cd 9 ae 9 M + 9 cd 9 af 9 be + 9 ac 9 be 9 df + 9 ac 9 M 9 de 



+ 9 bd 9 ae 9°^ + 9 bd 9°"^ 9° e + 9 ad 9 be 9°^ + 9 ad 9 b ^ 9° e + 9 b ° 9 ae 9^ + 9 b ° 9 a ^ 9 de 



(2.32) 



We just show the identity (2.30), the proof of the remaining ones is similarly. Let us 
consider the integral 



I{9cd) 



1 Tfl 

dqexp [ —Aq c g cd Aq L 



m 



27r i ehg 



(2.33) 



with the g c d as free parameters. Differentiation with respect to one of the parameters 
gives on the one hand: 



d 



dg a , 



I(9cd) 



im f 

2ehJ 



im 



dqexp [ ^-M c 9cdAq a )Aq a Aq b 



(2.34) 



and on the other 



dg, 



ab 



d * . / m d 1/9 1 / m 
—I(g cd ) = \ n _ : ^ a _ 9 1/2 = 



27r i eh dg ab 



+ 0-V2 ab = _ 1 abj, y / 2 go 

2 V 2irieh y y 2 y Ky 1 v ; 



Here we have used the formula for the differentiation of determinants: dg/dg a b = g g ab . 
Combining the last two equations yield (2.30). 

Let us denote by £ = q" — q' and q = q" '. The various Taylor expanded contributions 
yield: 



9 1/4 (q-09 1/2 (q-0^9' 6/ \q) 



3/4/ 



a/-b 



i-r a e a + -(4r a r 6 + 3r a , 6 )e a £ 



(2.36) 



exp 



i m 

2eh 



9ab(q-m a e 



exp 



im 

2eh 



ar-b 



9ab(q)t a t 



X 



1 + 2\eh^ ab ^ bcd ^ a ^ C ^ d ~ 8ieh if /av ^ c ' d dau^vdFbc + 9uv^ad^bcj C°C 6 C°C d 

(2.37) 



+ 9a V 9d u ri c T u ef Ce^ d CC f 



Here the various derivatives of the metric tensor g ab have been expressed by the 
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Christoffel symbols. Thus combining the last two equations yield: 



S 1/4 (9-0s 1/2 (4-6exp 



1 m 



2eh 



g ab ( q - oee 



g 3 '\q)es V 

2 



X 



1 / m 



i m 



9a b (q)d a e 



2 V 2 i eh 



2efr 

9av gdu vi c v u ef i a ei c i d i e i f 



- {gavTlcd + 9a U r a vd T% c + 9 u V r u ad T v bc )e^ d + l(^ a T b + 3T a , b )CZ 

(2.38) 

Let us start with the S a b-terms. We get immediately by equation (2.30): 



Similarly: 



B ab = - ^V. 
2rn 



B a =- 1 -^[T a g ab + (d a g ab )] 



For the £ 2 - and £ 4 -terms in So we get: 

^(4r a r 6 + 3T aib )ce~g ab (^ a r b + 3r a , 6 ), 

~~ g • ^ (dav^bc,d + dau^vd^bc + 9 uv^ad^bcj C° ^ 'C C 

h[ &T a T b + T c abc + 2r ajb 



(2.39) 



(2.40) 



(2.41) 



8m 1 



+ g U v9 cd (2K d rz c + r u ab r v cd ) + 5r ab r c bd + 2r d ac r 

For the £ 6 -terms equation (2.32) yields 

2 

g av g d urt c T u ef t a et% d t e t f 

~AT a T b + 4T ab T c + 4T d ac T c bd + g uv g cd (2T u ac Tl d + r u ab T v cd ) 



bd 



(2.42) 



If m 



2 V 2 i eh 



ieh ab 

~8m y 



(2.43) 



Therefore combining the relevant terms yields finally: 



lm 

2eh 



9ab(q-m a d 



J g 1/4 (q-£)gHq-0 exp 

. / m \~ D / 2 ( ie 
~ \2-K\eh) 



exp ( —AV Wey i 



(2.44) 
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Inserting all the contributions into equation (2.67) yields the Schrodinger equation 
(l-l)- 

Let us emphasize that the above procedure is nothing but a formal proof of the 
path integral. A rigorous proof must include at least two more ingredients 

1) One must show that in fact 

lim \[e- iTH/h -K(T)]^\ ^0 (2.45) 

N—toc 

for all \1/ G H (H: relevant Hilbert space). 

2) One must show that the domain T> of the infinitessimal generator of the kernel 
K(T) is in fact identical with the domain of the Hamiltonian corresponding to 
the Schrodinger equation (1.1), i.e. the infinitessimal generator is the (selfadjoint) 
Hamiltonian. 

It is quite obvious that these strong mathematical requirements hold only under certain 
assumptions on the potential involved in the Hamiltonian. Here Nelson [80] has shown 
the validity of the path integral for one dimensional path integrals. Some instructive 
proof can be furthermore found in the books of Simon [87] and Reed and Simon [85] . 
Also due to Albeverio et al. [1-3] is a wide range of discussion to formulate the Feynman 
path integral without the delay of go back to the definition of Wiener integrals. 



3. Product-Ordering 

In order to develop another useful lattice formulation for path integrals we consider 
again the generic case [43] . We assume that the metric tensor g a b is real and symmetric 
and has rank(g ab ) = D, i.e. we have no constraints on the coordinates. Thus one 
can always find a linear transformation C : q a = C a bVb such that Cci = ^A a by a y b 

with A a b = C ac g c dCdb and where A is diagonal. C has the form C a b = ui^ where 
the "Ur 3 ' (be {1, . . . , d}) are the eigenvectors of g ab and A a & = fl5 ac bb c where / o 2 ^0 
(a e {1, . . . , d}) are the eigenvalues of g a b- Without loss of generality we assume f% > 
for all a e {1, . . . , d}. (For a time like coordinate q a one might have e.g. f% < 0, but 
cases like this we want to exclude). Thus one can always find a representation for g a b 
which reads, 

9ab(q) = h ac (q)h bc (q). (3.1) 

Here the h a b = C ac f c C c b = ui" 1 f c Uc are real symmetric D x D matrices and satisfy 
h a bh bc = 5a- Because there exists the orthogonal transformation C equation (3.1) 
yields for the y-coordinate system (denoted by M y ): 

Kb(y) = fcWachc- (3.2) 

equation (3.2) includes the special case g a b = h-ab- The square-root of the determinant 
of 9abi \fg and the Christoffels T a read in the q-coordinate system (denoted by M q ): 
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The Laplace-Beltrami-operator expressed in the h ab reads on M q : 



Afj = I h ac h bc 



d' 2 



QqaQqb 



+ 



P)h ac P)h hc h 

un fobc _|_ if-c u ' b _|_ h ac h bc 
dq a dq a h 



_d_ 

dq l 



(3.4) 



and on M v 



J a 



+ 



j> ^Ja,a I 

h J dy a 



(3.5) 



With the help of the momentum operators (3.3) we rewrite the Hamiltonian in the 
" product-ordering" form 



H = "£ A ^ + V{q) = ^ haC ^P^h bc (q) + V(q) + AV Prod (q), 
with the well-defined quantum correction 



(3.6) 



AV, 



Prod 



8m 



Ah ac h bc ab + 2h ac h oc '^ + 2h 

,ao h 



acrbc h,ab 



ac I j^bc 



■b- 



be 



h 



' a h 



-h a °h 



ac-ibc h,ah } b 



h 2 



On M y the corresponding AVp rod (y) is given by 



(3.7) 



AV Prod (y) = 



8m P a 



fb,a 

X 



fa,c 



fa 



+ 4^ 



fa, a ( 2 fa, a fb,a 
fa V fa fb 



+ 2 



fb,a 

X 



. (3.8) 



The expressions (3.7) and (3.8) look somewhat circumstantial, so we display a special 
case and the connection to the quantum correction AVw ey i which corresponds to the 
Weyl-ordering prescription. 

1) Let us assume that A ab is proportional to the unit tensor, i.e. A ab = f 2 S a b- Then 
AVp ro d(y) simplifies into 



AV Prod {y) = h 



2 D-2(4-D)f% + 2f-f, aa 



8m 



f 4 



(3.9) 



This implies: Assume that the metric has or can be transformed into the special 
form A a }, = f 2 5 ao . If the dimension of the space is D = 2, then the quantum 
correction AVp rod vanishes. 
2) A comparison between (3.7) and (2.3) gives the connection with the quantum 
correction corresponding to the Weyl-ordering prescription: 



h 



AV Prod = AV Weyl + — (2h ac h bc , ab - h^ a h% - h™ b h\) . 



In the case of equation (3.2) this yields: 



(3.10) 



^ fa, a fafa, 



AV Prod (y) = AV Weyi + 



Am 



f 4 

J a 



(3.11) 
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These equations often simplify practical applications. 

Next we have to consider the short-time matrix element < q"\ e~ lTH \q' > in or- 
der to derive the path integral formulation corresponding to our ordering prescription 
(3.6). 

We proceed similar as in the previous section. We consider the short-time approxima- 
tion to the matrix element [e = T/N, g^' = g(q^')]: 



< qU)\ e - ieH / h \qV- iy > >~< q {j) \l - ietf/fc^-i) > 



D 



e ipAq ^' dp 



(2tvJi) 

- ^ < q (j) \h ac p aPb h bc \qti-V > -11 < q&\V + AV Prod \qU-V > . (3.12) 
2m n 



Therefore we get for the short-time matrix element (e < 1): 



< Q 



(i)| 



x exp 



itH/h I (j-i) > _ 



*frF ^ 2mh 



[9 U) 9 U ~ 1} 
(2nh) 



D 



dp 



h ac {q ( j))h bc {q u-l) )papb _ Ji V ( q (j)) _ j-AV Prod {q^) . 

(3.13) 



The choice of the "post point" q^' in the potential terms is not unique. A "pre- 
point", "mid-point" or a "product-form" -expansion is also legitimate. However, 
changing from one to another formulation does not alter the path integral, because 
differences in the potential terms are of 0(e), i.e. of 0(e 2 ) in the short-time Feynman 
kernel and therefore do not contribute. The Trotter formula e - lT ( A + B ) : = s _ 
\im.]\i^ 00 (e~ lTA / N e~ lTB / N ) N states that all approximations in equations (3.12) to 
(3.13) are valid in the limit N — > oo and we get for the Hamiltonian path integral 
in the "product form" -definition [h^ = h a b(q^)]: 



N-l . N , (j) 

K(q",q';T) = W) 9 {q")\^ Jim ]J / dq^ X ]J 

f • N r 

X exp < 1 V pAq^ - —h ac ^h bc ^- l) p^p[ j) - eV(q^) - eAV Prod (q^) 
I n 2m 



(3.14) 



Performing the momentum integrations we get for the Lagrangian path integral 
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in the "product form" -definition: (PF=Product-Form) 
K(q",q';T) 



«(*")=«" 



- I 



^fgV PF q{t) exp <j - 



rn 



-h ac h bc q a q b - V(q) - AV Prod (q) 



dt 



q{t>)=q> 



ND N-l , 

* n / \A^w<> 



N 



xexp|^ Aq a ^Aq b ^ - eV(q^) - eAV Prod {q^) 



(3.15) 



In the last step we have to check that the Schrddinger equation (1.1) can be deduced 
from the short-time kernel of equation (3.15). Because one can always transform from 
the q-coordinates to the y-coordinates, which is a linear orthogonal transformation 
and thus does not produce any quantum correction in the path integral (3.14) defined 
on M y , we shall use in the following the representation of equation (3.2). Looking 
at the Hamiltonian formulation of the path integral (3.14) we see that the unitary 
transformation y = C~ x q changes the metric term p a h ac h cb pb — > p a A- ah Pb: i.e. in the 
correct feature, the measure dq^dp^ — > dy^'dp^' due to the Jacobean J = 1 so 
that the Feynman kernel is transformed just in the right manner. 

We restrict ourselves to the proof that the short-time kernels of equations (2.25) 
and (3.15) are equivalent, i.e. we have to show (y = (y" + y')/2): 



[g(y')g(y")} ^v^M ex P 



1 777 1 f 

—A ab (y)Ay a Ay b - -V(y) -ie%AV w {y) 



exp 



l ^Uy')fa{y")A 2 y a - jV(y") - iehAV Prod (y")} \. (3.16) 



Clearly, e -'^ v ^)/ h = e -^v(v")/^ for the potential term. It suffices to show that a 
Taylor expansion of the g and the kinetic energy terms on the left-hand side of equation 
(3.15) yield an additional potential AV given by 



AV{y) = AV Prod (y) - AV Weyl (y) = 



fl 2 fa,a(y) ~ fa(y)fa,aa(y) 



4m 



ft{y) 



(3.17) 



We consider the (7-terms on the left-hand side of equation (3.16) and expand them in 
a Taylor-series around y' . This gives (£ a = (y" - y' a ), f a (y') = f a ): 



[g{y')g{y")\ ^v^I- 



i _ 1 fcfc,ab fc,afc,b C g C b 

8 P c * \ 



(3.18) 



Exploiting the path integral identity (2.30) we get by exponentiating the (3(e)-terms, 



[g{y')g{y")\ ^V^y-exp 



i eh fafa,bb fa,b 



8m 



f2 f2 
J a J b 



(3.19) 
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Repeating the same procedure for the exponential term gives: 



exp 



lm 

2eK 



exp 



lm 

2eh 



fa{y')f a (y")CC 



X 



i e 



1 (fcfc,ab fc,afc,b) £, £, £, £, 



(3.20) 



We use the identity (2.30) to get 



exp 



lm 

2eh 



a <-b 



exp 



i m 



2eh 
x exp 



fa(y')fa(y")ee 



lefl fafa,bb fa,b ^ \tfl fafa,aa fa, a 



8m 



f2 f2 
J a J b 



Am 



(3.21) 



Combining equations (3.20) and (3.21) yields the additional potential AT^ and 
equation (3.16) is proven. Thus we conclude that the path integral (3.15) is well- 
defined and is the correct path integral corresponding to the Schrodinger equation 
(1.1). 

A combination of lattice prescriptions were in the metric g a b in the kinetic term in 
the exponential and in the determinant g different lattices are used exist also in the 
literature. Let us note the result of McLaughlin and Schulman [70] : 



K W >, T) = w ta (^) Nm n / 

7 = 1 



x exp 



. N 

-Y 



h 



3 = 1 



m 
Ye 



g ab (q^)Aq^ a q^ b -eU(q 



0> 



with the effective potential 

U(q) = V(q)--A(q). q -F(q), 
c 

A(q) an additional vector potential term and F(q) given by 

771/ \ f-2 j-i I ab cd , ac bd , ad bc\ 

F(q) = -n F ahcd [g g +g g +g g J 

F a bcd i^Qyyi (y9ab,cd 2(7 J 'T abeX cdf^j 



(3.22) 



(3.23) 



(3.24) 



[r a bc = \{9ab,c + g ac ,b - 9bc,a)]- According to D'Olivio and Torres [21] the effective 
potential can be rewritten as 



U(q) = V(q) - -A(q) .q + — R + g ab (r c ad r d bc + d b F c ac ) 
c 8m L 

which is up the derivative term the result of Mizrahi [72] . 



(3.25) 
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4. Space-Time Transformation 

Let us consider a one-dimensional path integral 

x(t")=x" 



K(x",x';T) = 



Vx(t) exp 



x{t')=x' 



^x 2 -V(x) )di 



(4.1) 



It is now assumed that the potential V is so complicated that a direct evaluation of 
the path integral is not possible. We want to describe a method how a transformed 
path integral can be achieved to calculate K{T) or G{E), respectively, nevertheless. 
This method is called "space-time" transformation. This technique was originally 
developed by Duru and Kleinert [27,28]. It was further evolved by Steiner [90], Pak 
and Sokmen [82], Inomata [59], Kleinert [62, 64] and Grosche and Steiner [49]. In the 
rigorous lattice derivation of the following formulas we follow ourselves [49] and Pak 
and Sokmen [82]. 

Let us start by considering the Legendre-transformed of the general one-dimensional 
Hamiltonian: 

a 2 



which is hermitean with respect to the scalar product 

C/1,/2) = J h(x)mx)J(x)dx, J{x) = J Xv ^ dx ' 
Introducing the momentum operator 



hfd l rM1 
ft= iU + 2 r(l)1 ' 



F(x) = 



d\nJ(x) 

dx 



He can be rewritten as 
with the corresponding path integral 



2 



+ V(x) + — \r 2 (x) +2T'(x) 



E 



K E {x", x'- T) = e iTE / h K{x", x'- T), 

where K(T) denotes the path integral of equation (4.1). 
Let us consider the space-time transformation 

x = F(q), dt = f(x)ds 

with new coordinate q and new "time" s. Let G(q) = T[F(q)], then 



H E = 



h 2 1 

2m F' 2 {q) 



+ V[F(q)] - E. 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 
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2 

With the constraint f[F(q)] = F' (q) we get for the new Hamiltonian H = JHe'- 



H 



h 2 


' d 2 


2m 


dq 2 


h 2 


- d 2 


2m 


dq 2 



d 



F'(q) J dq 
+ f[F(q)]\V(F(q)) - E], 



+ f[F(q)W(F(q))-E] 



(4.9) 



where T(q) = G(q)F'(q) — F"(q)/F'(q). The corresponding measure in the scalar 
product and the hermitean momentum are 



h 



d 1 

dq + 2 



(4.10) 



The Hamiltonian H expressed in the position- q and momentum operator p q is 

6 = It + W^™)^ + AV W ( 4 - U ) 
with the well-defined quantum potential 



F"(q) 
F'(q) 



2 J ^ + (G(q)F>(q)) 2 + 2G'(q)F>(q) 



(4.12) 



Note that for G = 0, AV is proportional to the Schwartz derivative of the transfor- 
mation F. The path integral corresponding to the Hamiltonian H is 



l(t") = q' 



K(q'\ q'; a") 



J Vq{s) 



q(t>)=q' 

x exp 



i 

^ Jo 



m 



q 2 -f[F(q)}[V(F(q))-E]-AV(q) 



ds } . (4.13) 



As it is easily checked it is possible to derive from the short time kernel of equation 
(4.13) via the time evolution equation 



*(«";«) = / K(q",q'-s")*(q'-s)dq' 
the time dependent Schrodinger equation 



H^(q;s) = ih—^(q;s). 



(4.14) 



(4.15) 



The crucial point is now, of course, the rigorous lattice derivation of K(s") and the 
relation between K(s") and K(T). It turns out that K is given in terms of K(T) by 
the equations 



1 r°° 

K(x",x';T) = ^rj- J e- iTE / n G(x",x';E)dE 



(4.16) 



G(x",x';E) = i[f(x')f(x")}* I K(q",q';s")ds". 
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This we want to justify! 

Let us consider the path integral K(T) in its lattice definition 



K(x",x';T) = lim 



m 



n^oo \2irieh 



ND/2 



N-l 



n / 



N r 



x exp 



-y 



h 



3 = 1 



B( x U)- x U-i)f- eV ( 2 U)) 

2(: 



(4.17) 



We consider a D-dimensional path integral. To transform the coordinates x into 
the coordinates q we use the midpoint prescription expansion method. It reads that 
one has to expand any dynamical quantity in question F(q) which is defined on the 
points q(fi and q^' -1 ) of the j^-interval in the lattice version about the midpoints 
q(fi = \{q^ + q^~^) maintaining terms up to order (q^ — q^' -1 )) 3 . This gives: 

AF(q^) ee F{q&) - F(q^) = f(qW + ^) - f(qW - ^) 



dF(q) 



dq r 



+ -Aq^Aq^Aq 



U) d 3 F(q) 
dq m dq n dq k 



+ 



Introducing the abbreviations 

' m ~ dq m 



q=q 



U) 



7 U) 

,mnk 



d 3 F(q) 



dq m dq n dq k 



q=q 



U) 



we get therefore 

A 2 F^ = 



12 V m J 

according to equation (2.30) with the abbreviation 

F-LM U) ) = (F$Ftp)-\F^F^ 

+ (F^F^)-\F^F^ + (F$Ff)-\F^F$)-\ 

Furthermore we have to transform the measure. Because 



N-l 



N-l 



n dxU) = ii 

3 = 1 3 = 1 



dF(q 



dqti) 



dq 



N-l 
3 = 1 



(4.18) 



(4.19) 



~ Aq$Aq^F$F$ + 1 Aq£ Aq™ Aq^ Aq® F<£ F% 
=Aq^ Aq^F$F$ + 1 (^) ' F%F~Uq U) ) (4-20) 



(4.21) 



(4.22) 
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does not have a symmetric look, we must expand about We get: 



N-l 



n cfaw = [F, q w)F-,M')Y l/2 n d( i u) n to^)^- 1 )] 



N-l 



N 



N-l 



N 



3 = 1 



3 = 1 



X 



1 _ i _ F ;g>roB (gQ-)) \ m m 



N-l N 



x exp 



8m 1 ' m ' n } 



F; q (Q U) ) 



(4.23) 



Thus we have the coordinate transformed path integral 
K(q",q';T) 



3 = 1 



N , . 



3 = 1 



-AqV>AqWF im (qW)F in (qW) - eV(q^) 



_ ^( F (3) F (3h-l f F- ;q , m (Q U) )F; q AQ U) ) _ F; q ,mn(q^) 

8m l ' m j V FVqU)) Fq {qU)) 



eh 2 
/. 

8m ' m > n ^' mnkl 



F U) F U) F -i ( 



(4.24) 



This path integral has the canonical form (2.24). In order to see this note that 
the factor F m (gW')F n (^) in the dynamical term can be interpreted as a metric 
9nm appearing in an effective Lagrangian. Thus we can rewrite the transformed 
Lagrangian in terms of this metric (Gervais and Jevicki [38]): Firstly, we have 
F-, q (q) = a/ det(g nm (q)) = \/g(q). In the expansion of the determinant about the 
midpoints we then get 



[f-M^f-M'- 1 )] 1 ' 2 

1 



q{q^){l + 



16 



(4.25) 



Using now successively the identities 

' (4.26) 
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we can cast the quantum potential into the form 



AV(q) = ^ mn (q)rL(Qn m (Q) 



(4.27) 



which is just the Weyl-ordering quantum potential without curvature term. In the 
present case the curvature vanishes since we started in flat Euclidean space which 
remains flat during the transformation. 

In particular equation (4.24) takes in the one-dimensional case the form [F' = 
dF/dq}: 



K( q '\ q ';T) = [F>( q ')F'( q ")}--> Jim 



m \ N / 2 



N-l 



N^oo \27rieh 



n / 



N 



m 

Mi 



Aq$Aq^F' 2 (q^)-eV(q^) 



eh 2 F' 2 {q^) 
8m F' A ( q ^)) 



(4.28) 



which is the same result as taking for D = 1 in AVweyl the metric tensor g a i, = F' . 
The infinitesimal generator of the corresponding time-evolution operator is He- To 
proceed further we perform a time transformation according to 



dt = f(s)ds, s' = s(0 = 0, s" = s(t") 



(4.29) 



and we make the choice f = F' . Translating this transformation into the discrete 
notation requires special care. In accordance with the midpoint prescription we must 
first symmetrize over the interval — 1) to prefer neither of the endpoints over the 
other, i.e. 



Expanding about midpoints yields 



€ ~8^F >2 ( q ^) \l + 



F"'{qW) _ ( F"( q ^) \ 2 



(4.30) 



(4.31) 



Insertion yields for each j^-term [identify V(x^) -> V[F(qW)] = V(q^)] 



N N-l 

m \ t 



N 



2n i eh 

N-l 



exp < - 



h 



rn 
2~e 



AV j) -eV(x {j) ) + eE 



N-l N 

nnwn(^) 



i f m 

l .77\27 



F ,2 (qW)A 2 qW + ±-F'( q ^)F" , ( q ^)A 4 q {j) 

X. £ 



eV(x U) )+eE 
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N-l 



N 



3=1 

i | m 



m 



TV 

x ]^[exp 



x 



ft 1 250') 



| A 2 gO) ( F"'(q^) ( F"(qW) } 

i v /- 7 (//' ; 



12 F'(gO)) 

4 iV-l 



(qO)) \F'(q(i)) J 
d^F' 2 {q^)\v{q i3) )-E 



F'{q')F'{q") 



m 



2ni5h 



i=i 



x exp 



At I ZiO L 



8m 



F'(gO)) y F'(gO)) 



(4.32) 



Let us assume that the constraint 



f dsf[F(q(s))] = T 
Jo 



(4.33) 



(4.34) 



has for all admissible path a unique solution s" > 0. Of course, since T is fixed, the 
"time" s" will be path dependent. To incorporate the constraint let us consider the 
Green function G(E): 

poo 

G{x", x'; E)=i dT e iTE / h K(x", x'; T) 
Jo 

K(x",x';T) = ^-j- J dE e- iTE/h G(x",x';E). 

This incorporation is far from being trivial, in particular the fact that the new time- 
slicing <j0') is treated in the same way as the old time-slicing e, i.e. £0) = S is considered 
as being fixed, but can however be justified in a more comprehensive way, see e.g. [6, 
35, 59]. We now observe that G(E) can be written e.g. in the following ways [62]: 

/•OO 

G{x",x';E) = f(x") / <x"\e- ia " f ^ H - E ^ h \x' > ds" (4.35a) 

Jo 

/•OO 

= f(x') / <x"\e- is "^ H - E ^/ h \x' >ds" (4.35b) 
Jo 



y/f(x')f(x") ^ < x »\ e -""yffWW-E)yfJG)/h | x , > ds n 
Jo I 



(4.35c) 



The last line give for the short-time matrix element in the usual manner 

<x ^)\ e -^ {3) f^{x)(H-E)n{ x )/n ,(i-l) > 
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=< x^\x 



(/)|r0 _ 1) > . v ' /(,•:./:) /,,.,/)) < X U)\H-E\xV-V > 



< F(q^)\F( q U-V) > - 1 Jf[F(qU))]f[F(qU))} < F(q^)\H E \F(q 



h 



> 



iS 



< q^lqU- 1 ) >--r< q {3) \H\q^~^ > 



2tvJ 



h 



dp„ exp 



m 



exp 



i m 



2Sh A 2 q (i) ~ jrF ,2 (jW)(Vtf») -E)- 5 -AV{q^) 



2n\5h 

where the well-defined quantum potential AV is given by 



AV(q) 



1_ 

8m 



z F»'(q) jF"(q) 



F'{q) 



F'{q) 



(4.36) 



(4.37) 



Note that this quantum potential has the form of a Schwartz-derivative. This 
procedure justifies our symmetrization prescription of equations (4.30). Thus we arrive 
finally at the transformation formulae: 



i r°° 

K(x",x';T) = — — / dEe- iTE / h G[F(q"),F(q'y,E} 
Zirih J_ 00 

poo 

G[F{q"),F{q')-E] = \[F'{q")F\q')}-^ / ds" K(q" ,q>; s") 

Jo 

where the path integral K(s") is given by 

N-l f 
II / dq<» 



(4.38) 



K(q ,q;s ) = lim 



iV^oo \ 27T i 8Ti 

N r 

xexpU^ S A V j) - SF' 2 (q^)(V(q^) - E) - SAV(q^) 

j = l <- 



(4.39) 



This is exactly the path integral (4.13) and thus we have proven it. Note that only 
for D = 1 we have the simple transformation property 



2e 
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(4.40) 



For higher dimensional problem this may only possible for one coordinate kinetic term. 
Anyway, the requirement for higher dimensional problems are 



F = F F 



(4.41) 
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for at least one pair (m, n) £ (1, . . . , D). 

Finally we consider a pure time transformation in a path integral. We consider 
equation (4.35c), i.e. 



G(x", x'; E) = y/f(x')f(x") / ds" < q"\e 

Jo 



i8"y/f(H-E)y/f/h 



(4.42) 



where we assume that the Hamiltonian H is product ordered. Then 
G(x", x'\ E) 



y/ fix') f{x") r ds" H ( dq& f[ < q 

J ° 3 = 1 J 3=1 

f ds" lim FT / dq (j 



U)\ e -iSy/f(H-E)y/f/h\ q (j-l) > 



_ y/ fix') fjx" 



N - • ,0") 



x 



II / 7^d ex P ( i\ P U) *Q U) - 8y/f(gW)fW-») 



x 



h ac (qti-V)h bc iq^) 



p U) p U) + v ^ q U)) + AV Prod (q^) - E 



2m 



^') i(1 - D/ \-(^) ND/2 n/^ 



(J) 

2e VIW^WW) 



iqU)) 



-iefiqW)(yiqW) + AV Prod iqW)-E 

poo 

= {rn Ui-D/2) / K{q",q'-s")ds" 
Jo 



(4.43) 



with the path integral 



fj.ff\ if 

K(?",<zV) = / V$V PF q(t) 

Q(t')=q' 



x exp < — 

ft Jo 



-h ac h cb q a q b ~ f ( + AV Prod (q) - E 



ds } . (4.44) 



Here denote h ac = h ac /V7 and \fg — det(/i ac ) and equation (4.44) is of the canonical 
product form. Note that for D = 2 the prefactor is "one" . 
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5. Separation of Variables 

Let us assume that the potential problem V(x) has an exact solution according to 
J Vx(t)exp lj (^&-V{xY)dt = JdE x e- iE ^ n n(^x(x")- 

x(t')=x' L J 

(5.1) 

Here f dE\ denotes a Lebesque-Stieltjes integral to include discrete as well as 
continuous states. Now we consider the path integral 



K(z", z', x", x'; T) 
*<*">=*" d 



x(t")=x" 



I f d (z)f[g i (z)Vz i (t) J f[Vx k (t) 

f\ / i — 1 „ i j.1 \ / /c = l 



z(t')=z' 



x(t')=x 









m I 




7 ' 



d' 



i=l 



ND N-l 

m \ — 



k=i 



V{x) 
IK*) 



+ W(z) + AW(z) 



dt 



N^oo \ 2-K\eh 

. N 



= Hm (— ) ■ n //wn^w/n^ 

j = l ^ i=l ^ fc=l 



x 



m 
2^ 



i = l 



k=l 



— € 



V(x 



p( z U) 



^L + w{z u) ) + AW{z u)^ 



(5.2) 



Here (z, x) = (zi, x k ) (i = 1, • • • , d'; k = 1, . . . , d, d' + d = D) denote a D-dimensional 
coordinate system, gi and / the corresponding metric terms, and AW the quantum 
potential. For simplicity I assume that the metric tensor g a jj involved has only 
diagonal elements, i.e. g ab = dia,g(g j (z), g%(z), . . . ,g#(z), f 2 (z), . . . , f 2 (z)]. Of course, 

det(<7 a f,) = f 2d rii=i 9i = f 2d( ~*( z )- The indices i and will be omitted in the following. 
We perform the time transformation 



s = 



da 



s" = s(t"), s(t') = 0, 



(5.3) 



where the lattice interpretation reads e/[f(z^~ 1 ' ) )f(z^' > )] = 5^ = 5. Of course, we 
identify z(t) = z[s(t)] and x(t) = x[s(t)]. The transformation formulas for a pure time 
transformation are now 



1 f°° 

K{z", z', x", x'-T) = — — / dE e" 1 ET/h G(z" , z', x", x'; £) 



(5.4) 
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G(z",z',x",x';E) = i[f(z')f(z")} 



f/\]l-D/2 



ds"K{z",z',x",x';s"), 

(5.5) 



where the transformed path integral K(s") is given by 
K(z", z' , x", x'; s") 

( n\ " t H\ I' 

z(s )—Z X{S )=X 

z(0)=z' x(0)=x' 

m f g 2 (z) 



x exp 



, v f2{ y) z 2 + J - V{x) - f 2 (z) (W(z) + Ali o) + f{z)E 
J dE x e- [E ^"/ h ^l(x')^ x (x")k(z",z';s"), 



ds 
(5.6) 



with the remaining path integration 

z(s )=z 

K(z",z';s")= ^§^Vz{s) 



z(0)=z' 



f 2 ~ d (z) 



Xexp ^y 



m g 2 (z) . 



z 2 - f 2 (z) (W(z) + AW(zj) + f 2 (z)E 



ds } . (5.7) 



Of course, in the path integrals (5.6,5.7) the same lattice formulation is assumed as in 
the path integral (5.2). Note the difference in comparison with a combined space-time 
transformation where a factor {f{z')f(z")] 1 ^ would instead appear. We also see that 
for D = 2 the prefactor is identically "one" . We perform a second time transformation 
in K(s") effectively reversing the first: 



a= f f[z(u)]du, a" = s" 
Jo 



(5.8) 



with the transformation on the lattice interpreted as = 8^f(z^ 1 ')f{z^). 
Therefore we obtain the transformation formula? 



1 f°° 

K(z",z';s") = -^— / dE' e- iE ' a "/ h G(z",z';E') 

POO _ 

G(z",z';E>) = i[f(z')f(z")] E ^- 1 / da"e iE °"/ h K(z",z';a") 

Jo 

with the transformed path integral given by 

z(a )=z 

k{z",z';a")= J ^G\z)Vz(a) 



(5.9) 
(5.10) 



z(0)=z' 



x exp 



ft Jo 



rn 



g 2 (z)z 2 -W(z) - AW(z) + 



E' 

1W) 



da 



(5.11) 
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Plugging all the relevant formulae into equation (5.5) yields 



K(z",z>,x",x>;T) = \f(z')f(z")]- D/2 J dE x 9l(x')9 x (x") 



x 



1 f°° ,i r 

J J_ 



2nh 

z(o")=z 



dEe - iE (a"-T)/n 



2nh 



dE' / ds" e- i8 "^+ E '^ h 
Jo 



x 



/ ^ 



z) Vz(a) exp 



z(0)=z' 



^g 2 (z)z 2 -W(z)-AW(z) + ^ 



da 



(5.12) 



The <ia // <iE'-integration produces just a" = T, whereas the (iE'ds^-integration can be 
evaluated by giving E\ + E a small negative imaginary part and applying the residuum 
theorem yielding E\ = —E'. Therefore we arrive finally at the identity [45] 



K(z", z', x", x'; T) = [f(z')f(z")]- D / 2 J dE x n(x')*x(x") 



z(t")=z" 



x 



J y/G{z)Vz(t)e^r-^ 



z(t')=z' 



^g\z)z 2 -W(z)-AW(z)-^- ) 



dt > . 



(5.13) 



Note that this result can be short handed interpreted by inserting 

D/2 



m 



exp 



i / m 



h V 25^') 



T AV J) - 5 {3) V{x {3) ) 



= J dE xU) e-^xO)* ^^*^^-!))^^)) (5J4) 

with 5^^ = e /[f (z^ f (z^^)] for all j th and applying the orthonormality of the 
in each j t/l -path integration. Equation (5.14) describes therefore a "short-cut" to 
establish equation (5.12) instead of performing a time transformation back and forth. 
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III Important Examples 

In this Chapter I present some of the most important path integral solutions which 
are 

1) The free particle as a simple example. 

2) The harmonic oscillator in its basic form, where we allow in addition a time 
dependent frequency. 

3) Path integration in polar coordinates. We shall discuss the various features of 
properly defined path integrals including the "Besselian functional measure". 
In particular the (time-dependent) radial harmonic oscillator will be exactly 
evaluated. 

4) The Coulomb potential. 



1. The Free Particle 



For warming up we calculate the path integral for the free particle in an .D-dimensional 
Euclidean space. From the representation 

K(x",x';T) 

x(t")=x" 



x(t')=x' 



= lim 



m 



n^oo V27rie7i 



ND N—l „oo 



n 

j=l ->~°° 



exp 



N 



(1.1) 



it is obvious that the various integrations separate into a D-dimensional product. All 
integrals are Gaussian. Starting with (/i = 1, . . . , D): 

nOO r 

= f TT^r) / dx^ exp - -^-(4 2) - *W - ^-(4 1} - 4 0) ) 2 
\2niehJ J_ 00 F 2ieft v ^ ; 2ie% K ' 



m 



An\ Tie 



exp 



m 



(4 2) -4 0) ) 2 



(1.2) 



it is easy to show that 



m 



2tt i TieN 



exp 



m 



Tin V (a: " 



(TV) 



4 0) ) 2 



(1.3) 
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Thus we get in the limit N — > oo (note eN = T): 

D/2 



K(x",x';T) = 



m 



2niHT / 

2 r°° 



exp 



2ftT 



(x — X ) 



(27T) 



D 



exp 



ip(x"-a;') -iT 



ftp 2 
2m 



(1.4) 



dp. 



From this representation the normalized wave-functions and the energy spectrum can 
be read off (p G R D ): 



which is the known result. 

The corresponding energy-dependent Green function is given by (D = 1): 



(1.5) 



(x", x'- E )=iJ Q °° K (x", x'\ T ) ^ ET/h dT = y^^exp 



\x" — x'l 



h 



For the D-dimensional case we obtain 
G( D \x",x';E) 

= 2 i I I ! \.r - .r 



V2mEj . 
(1.6) 



2yrin 



2£' 



i(l-D/2) 



I \ 



h 



v / 2mE ) 



-(i) D/ W'--^) 4,1 "° /2 -i%(^^) 

where use has been made of the integral representation [40, p. 340] 



'(1.7) 
(1.8) 



dt t v exp ( — — — pt 



a 



= 2(- 

V4p 



K u {^fap), 



(1.9) 



and #„(*) = if e i7r ^ 2 M 1} (i^), tf±i(z) = v^T^ e " 2 • 



2. The Harmonic Oscillator 

After this easy task we proceed to the path integral calculation of the harmonic 
oscillator. We consider the simple one-dimensional case with the Lagrangian 

C[x, x] = —x 2 - -^x 2 + b(t)xx - e(t)x (2.1) 
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Here we assume that the various coefficients may be time-dependent. The path integral 
has the form 



x(t")=x" x(t")=x" 

K(x",x';t",t')= J Vx{t)e^ = J Vx(t)e*f* C[x ' ±]dt (2.2) 



x(t')=x' 



x{t')=x' 



In this quadratic Lagrangian, of course, an ordering problem appears. This is due to 
the stochastic nature of the path integral, the "time" -integral in the exponent is not 
a Riemann-integral but a so called "Ito-integral" . This we have already discussed in 
chapter II. The ordering ambiguity appears in the "b(t)"-term, where we have in the 
corresponding Hamiltonian a px-term. The consequence is that we have to use the 
mid-point formulation, i.e. (x^ = \{x^ +a;( J_1 )), and = c(t^), etc., see also 
Schulman's book [86]): 



K(x",x';t",t') 

o(t")-- 



Vx(t) exp 



x(t')=x' 



N 



xexp 



AT N—l „oo 



r ffh .9 CI b J O t/ \ • / \ 

— x — x + b{t)xx — e(t)x 



dt 



N r 



i v ^ m / 



3 = 1 



rj£> ( 3 ) X 



'') - e C —{x^f + b^x^(x^-x^) -v»'.i 



(2.3) 



In the following we use the notations in equations (2.2) and (2.3) as synonymous. Let 
us expand the "path" x(t) about the classical path xci(t), i.e.: 

x(t) = x C i(t) + y(t), 



where y(t) denotes a fluctuating path about the classical one. The classical path 
obeys, of course, the Euler Lagrange equations: 



SC\ 



xci,xci\ 



d dC[xci,x C i] dC[x C i,xci 



5x C i dt dx C i dx C i 

Expanding we obtain for the action 



= 0, x C i(t') = x', x C i(t") = x" cl . 

(2.4) 



S[x] = S[x C i] + - 



1 S 2 S\x] 



2 Sx 2 



x=xci 



y 2 = S[x C i(0,xci(t')}+ f 

Jt' 



1 1 v Q ° 2 



dt. 
(2.5) 
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The linear functional variation vanishes due to the Euler-Lagrange equations (2.4). 
For the path integral we now find 



K(x"x'; t", = exp ^S[x C i(t"),x C i(t')} jF(f, t'), 

y(t")=o t „ 
F(t",t')= J Vy(t)e,p ^ (f - u\t)y*) 



dt 



(2.6) 



Here we have used the abbreviations muj 2 (t) = c(t) + b(t), y' = y(t') and y" = y(t"). 
Now 



F(t",t') = lim F N 



= lim 

N^oo 



N 



x exp J2 [(y {j) - ^ J_1) ) 2 - ^ U)2 v U)2 



(2.7) 



Let us introduce aiV- 1-dimensional vector z = (x\, . . . ,xn-i) T and the (N — 1) x 
(N — 1) matrix B: 



/2-eV 1 ) 2 -1 

-1 2-eV 2 ) 2 



B = 





V o 


















\ 



2 _ e 2 w (iv-2)2 _ x 

-1 2-eV N - 1 ) 2 / 



(2.8) 



Thus we get 



iV 



AT 

m \ — 



2tt 



Therefore 

F(t",t') = 



m 



m 



m 



i 

\ 2 



2yrien det S 



27Tift/(tV) 

Our final task is to determine detS. Let us consider the j x j matrix 



/2-eV 1 ) 2 -1 

-1 2-eV 2 ) 2 





V o 













\ 





2-eV*'- 1 ) 2 -1 

-1 2-e 2 ^^ 2 / 



(2.9) 



where f(t",t')= limedetS. (2.10) 



(2.11) 
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One can show that the following recursion relations holds: 

det B u+l) = (2 - eV J ' +1) 2 ) det B^ - det B u ~^ 

with det = 2 - e 2 u^ 2 and det B^ = 1. Let us define ^ = edetB^, then we 
have 



,0+1) 



jl „(i-i) - _^2,,,(j + l)2 n 0-) 



Turning to a continuous notation we find for the function g(t) = f(t, t') a differential 
equation: 

g(t)+u 2 (t)g(t) = 0, with g{t') = 0, g(t') = 1. (2.12) 

The two last equation follow from <7(£') = go = lim e ^o edetS*- -* and g(t') = lim e ^o 
[g(t' + e) - g(t')]/e = lim e ^ (det - det B^) = 1. Finally we have to insert 
g (t") = f{t",t') into equation (2.10). 
At once we recover the free particle with g = t" — t' = T. 

The case of the usual harmonic oscillator with u(t) = to (time-independent) is given 
by ( 

g(t) = — sinu/T. 

To get the path integral solution for the harmonic oscillator we must calculate its 
classical action. It is a straightforward calculation to show that it is given by: 



S[xci(t"),x C i(t')} 



muj 



2 sin uT 

and we have for the Feynman kernel: 



(x'ci + x" 2 cl ) cosuT - 2x' cl x'ci 



(2.13) 



K(x",x';T) 



ex P i - 



27r i h sin uT J { 2ih 

By the use of the Mehler-formula [31, Vol.III, p.272]: 



(x' 2 + x' 2 ) cotuT — 2—. 



sincuT 



(2.14) 



n=0 



1 



exp 



Axyz- (x 2 + y 2 )(l + z 2 ) 
2(1 -z 2 ) 



(2-15) 

where H n denote the Her mite-polynomials, we can expand the Feynman kernel 
according to (identify x = ^muj/Ux 1 ', y = ^muj/Ttx" and z = e~ lujT ): 



K(x",x';T) = J2e-' lTE "/ h q* n (x')q n (x") 



n=0 



with energy-spectrum and wave-functions: 
E n = hu(n+ §), 



muj 



x exp 



muj 
~2~E 



x 



(2.16) 



(2.17) 
(2.18) 
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Equation (2.14) is as it stands only valid for < uT < it. Let us investigate it for 
larger times. Let 



nn 



T = — + t, with n G N ; < r < n/u (2.19) 
then we have sincuT = e 17rn sino;r, coscuT = e 1 7rn cos lot and equation (2.14) becomes 

K(x",x';T) = 



2-kU sin cur 



in f 1 \ imu 

x exp < + n H — — 

1 ] 2 \2 / 2hsmuT 



(x" + x' 2 ) cos ivT -2x'x" \ (2.20) 



which is the formula for the propagator with the Maslov correction (note sinter 
| sinu;T|). Now let r — > 0, i.e. we consider the propagator at caustics. We obtain 



K[x",x';—) - lim 



lo J t^o \ 2irihT 



m 



inn im 
° XP<5 + 2^ 



0r /2 + ; r /2 ) ^e-^xV 



lim 



exp 



-^-(-i)vr-i^ 



exp ^- ^m^x'- (-l)V') 



(2.21) 



Let us finally determine the energy dependent Green function. We use the integral 
representation [40, p. 729], a\ > a 2 , 9R(^ + n — v) > 0): 



coth 2iy — exp 
2 



x 



ai + a 2 



t cosh x 



-^2/x(^v /a i a 2 smhx)dx 



t^/aia 2 T(l + 2jj) 



W^( ai t)M^(a 2 t). (2.22) 



Here W Vjfl (z) and M v ^(z) denote Whittaker- functions. We reexpress e x = y/nx/2 
[Ii (x) +I_i(x)] apply equation (2.22) for v = E/2fko and use some relations between 
Whittaker- and parabolic cylinder-functions to obtain 



1 m 



G(x", x'; E) = — ,, . . 
v ' ' ; 2\ irhu \2 huj 



xD i 



2moo 



h 



(x' + x" + \x" -x'\) 



D i 



2moo 



(x' + x" - \x" -x'\) 



(2.23) 



From this representation we can recover by means of the expansion for the T-function 



r(«) = £ 



n=0 



-I) 



n\(z + n) J 1 



38 



III.2 



The Harmonic Oscillator 



and some relations of the parabolic cylinder-functions and the Hermite polynomials, 
the wave- functions of equation (2.18). 

Let us note that we can use equations (2.14, 2.23) to obtain recursion relations for 
the Feynman kernel and the Green function, respectively, of the harmonic oscillator 

in D dimensions. Let x G R n and define fj, = x' + x" , v = x' • x", then 



K( D \x",x';T) = ^^-K( D - 2 Hx",x';T) 



1 

= 27 6 



iu>T 



moo 



^ 2^jK^ D - 2 \x\x';T). (2.24) 



For the Green functions this gives 



G {D \x",x';E) = & D \ii\u- E) 
& D \^E) = ^& D ~ 2 \^E) 



(2.25) 
(2.26) 

(2.27) 



Introducing furthermore £ = \{\x' + x"\ + \x" — x'\), rj = \{\x' + x' n 



obtain for D = 1, 3, 5 ... . 



G^ D \x",x';E) = 
1 



1 m 



1 E \ ( 1 



r -- 



- 



V 2 ~ i 2 V ®i 



2 V Khu \2 hu J \2ir 
d 



D-l 

2 



-D-I + JL 



respectively, 



x 



moo 
~~h~ 



D-l 
2 



D D 



D-l 
2 



2mu 



V ' ' 1 2\j Txhoo \2 hu)\2n) 



D-l 
2 



2mu 



\x" — x'\) we 



2mu 



h 



v 



(2.28) 



E 

2 1 Su 



2muj 



(2.29) 



Let us note that the most general solution for the general quadratic Lagrangian is due 
to Grosjean and Goovaerts [52, 53]. 
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3. The Radial Path Integral 



3.1. The General Radial Path Integral 

Radial path integrals have been first discussed by Edwards and Gulyaev [30] and 
Arthurs [4]. Edwards and Gulyaev discussed the two- and three-dimensional cases, 
Arthurs concentrated on D = 2. Further progress have been made by Peak and 
Inomata [84] who calculated the path integral for the radial harmonic oscillator 
including some simple applications. See also [49] for the ordering problematics 
which give rise to the various quantum potentials in the formulation of the radial 
path integral. In this section we derive the path integral for D-dimensional polar 
coordinates. We follow in our line of reasoning references [49] and [92], where these 
features have been first discussed in their full detail. Similar topics can be also found 
in Bohm and Junker [5] from a group theoretic approach. We will get an expansion 
in the angular momentum /, where the angle dependent part can be integrated out 
and a radial dependent part is left over: the radial path integral. We discuss some 
properties of the radial path integral and show that it possible to get from the short 
time kernel the radial Schrodinger equation. 

The path integral in D-dimensions has the form: 



Now let V(x) = V(|x|) and introduce D-dimensional polar coordinates [31, Vol.11, 
Chapter IX]: 



x\ = r cos 6*i 

X2 = r sin6>i cos 62 

£3 = r sin 6>i sin 9 2 cos 6> 3 

xd-i = rsin6>i sin 62 ■ ■ .sm6o-2 coscj) 
xd = r sin 6*i sin 2 ■ ■ ■ s\ti6d-2 sin <f> 



K {D \x",x'-T) 




x(t')=x' 




(3.1) 
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where < 9 V < ir {y = 1, . . . , D - 2), < < 2tt, r = ^2u=i x l^ 
V(x) = V(r). We have to use the addition theorem: 

cos^ (1 ' 2) = cos^cosflf 



1/2 



Therefore 



D-2 



+ J2 «* e$ +1 cos 9% +1 J] sin 0« sin 9^ + f[ sin sin 9^ , (3.3) 



)(2) 



D-l 



m=l 



n=l 



n=l 



where ^A 1 ' 2 ) is the angle between two .D-dimensional vectors and so that 
X (X) . x ( 2 ) = r (!) r (2) Qosifji 1 ' 2 ), The metric tensor in polar coordinates is 



(g ab ) = diag(l, r 2 , r 2 sin 2 1? . . . , r 2 sin 2 ^ . . . sin 2 9 D - 2 ). 
If D = 3 equation (3.3) reduces to: 

cos-i/> (1 ' 2) = cos# (1) cos6> (2) + sin^ 1 ) sin#( 2 ) cos((j)^ — 4>^). 
The D-dimensional measure dx expressed in polar coordinates is 

D-l 



(3.4) 



dx = r^drdSl = r D ~ x \\ (sin 9 k ) D - 1 - k drd9 k 

k=i 

D-l 

dQ= Yl (sin 9 k ) D - 1 ~ k d9 k . 



k=l 



(3.5) 



d£l denotes the (D — 1) -dimensional surface element on the unit sphere S D 1 and 
Q(D) = 2n D / 2 /T(D/2) is the volume of the D-dimensional unit (S ' 1 -) sphere. The 
determinant of the metric tensor is given by 



g = det(g ab ) = (r " 1 J] (sin 9 k ) D - 1 ~ k ) 
^ fe=i ' 



(3.6) 



With all this the path integral equation (3.1) yields: 
K( D \r",{9"}y,{0'};T) = 

/ Tn \ N D /2 r°° C roo 



dQ 



(N-l) 



X 



II ex Pi ^0') + r 0-D -2r y) r (3 -_ 1) cos^- 1 )] - ^( 



j'=i 



{9} denotes the set of the angular variables. For further calculations we need the 
formula [40, p.980]: 



3 z COS Ip 



(3.8) 



1=0 
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(for some v ^ 0, —1, —2, . . . ), where C\ are Gegenbauer polynomials and 7 M modified 
Bessel functions. Equation (3.8) is a generalization of the well known expansion in 

three dimension where v = \ (remember Cj 2 = Pi, [40, p.980]): 



^Z COS 1p 



I oo 

V 2z 1=0 



(3.9) 



Note: It is in some sense possible to include the case D = 2, i.e. v = if one uses 
/lim A ^or(A)C z A = 2cos/V' [40, p.1030], yielding finally [40, p.970]: 



^Z COS 1p 



E J *(*) eI 



kip 



(3.10) 



k= — oo 



The addition theorem for the real surface (or hyperspherical) harmonics S^ 1 on the 
S^-sphere has the form [31, Vol.11, Chapter IX]: 



M 



^S^mrtV) = n 1 {D) 2l+ D D 2 2 C^(cos^) (3.11) 

with Q = x/r unit vector in R d and the M linearly independent Sj 1 of degree / with 
M = (21 + D - 2)(l + D - Z)\/l\(D - 3)!. The orthonormality relation is 



/ 



(3.12) 



Combining equations (3.8) and (3.11) we get the expansion formula 



e *(n (1) 'fi (a ») =e2 cos^( 1 - 2 ) = 2n 



2tt N 2 

z 



E E s?(nU)s?(nW)i l+£ ^(z). (3.i3) 

1=0 fi=l 



Let v = (D — 2)/2 > in equation (3.8), then the jth term in equation (3.7) becomes 



R J =ex V) 2&( r U) +r U-i))-J V ( r U))\ ex P 2i^ r (7) r 0-i) cos V' UJ J 



D-2 

2ieft \ 2 /D-2 



mr U) r U-i) 

oo 

X 

z,=o 



exp 



im / 2 2 \ ie Tr/ x 

2^( r «) + r (i-i))-n^-)) 



D 



/ 2yrien 
= 2?r - 



\ mr (j) r U-i) 



D-2 
2 



exp 
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With the help of equations (3.11) and (3.14) now (3.7) becomes: 
K( D \x",x';T) 



oo M 



^T^EE 5 f( ')5f(n") 

1=0 n=l 

N eQO 



x lim 



m 



N 



x y\ ex p 



i m 



2eh 



I* oo 

r (l) dr (l) / ^(JV-l)^(iV-l) 

Jo 

I l+ E_l[^r {j) r u _ 1) \. (3.15) 



Therefore we can separate the radial part of the path integral: 



K^ D \r", {$"}, r', {$'}; T) = Sl~\D) £ l ^-J. cf^ (cos / "))K,(r", r'; T) 

(3.16) 



with 

Ki D \r",r';T) = (r'r")" 



p-2 ( m \ 

2 lim - . - 



N—>oo 



n/2 r°° r°° 
/ dr^y / dr 
Jo Jo 



(N-l) 



X 



and the functional measure 

.(D) 



N [ . JV 

n M^ro-i)] • exp ^ ij] 



m 
27 



( r (j) - r (j-i)) 2 " el/ ( r (j)) 



(3.17) 



W (^(j)) = V 2nz (j) e 2(j)/ z+^(^))> 



(3.18) 



where = (m/ i e^r^r^^y 

In the literature often use is been made of the asymptotic form of the modified 
Bessel functions 



I„(z) ~(2nz)-^ e z ~^ 2 -^ /2z (\z\ > 1, > 0). 

Then the functional measure becomes (we ignore 9?(z) = 0): 



^ (D) ( 2 0')) - ex P 



i efo 



2mr 0) r _ 1) 



1 + 



D — 2 



(3.19) 



(3.20) 



and reads: 



K[ D \r",r'-T) = {r'r")- ki± lim 



N^oo \27r i eft 



N/2 



/•OO />C 

/ dr (1 y 
Jo Jo 



drtN- 



(N-l) 



X 



i ^ 



i=i 



^ 2 / 4 _ T7" / \ 

2e v ' w '^"^ - 2mr (j) r (j _ 1) e ^ r C0J 



. (3.21) 
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This last equation seems to suggests a Lagrangian formulation of the radial path 
integral: 



r(t")=r" 



k\ D \r",r';T)-- 



Vr(t) exp 



r(t')=r' 



2 2mr 2 



- V(r) 



(3.22) 

with k\ D \r",r';T) = ( r V')(^-i)/2^( D )( r V' ; T). But nevertheless, equation (3.17) 
can be written in a similar manner with a nontrivial functional measure: 



r(t")=r' 



k 



(D) 



(r",r';T)= J ^[r 2 ]Pr(t)exp^ 



r(t')=r' 



jf 2 - V(r) 



dt 



(3.23) 



and the whole /-dependence is in as defined in equation (3.18) in the lattice 
formulation. The asymptotic expansion of I v in equation (3.19) in problematic because 
it is only valid for $l(z) > 0, but we have $l(z) = 0. The complete expansions in this 
case reads [40, p. 962] 



(-l) k T{v + k+\) 



2tt*^ ( 2z ) k k\T(v-k+\) 



+ exp[-z± + g 1 T{v + k+\) 



V2 



V2 



nz 



exp 



nz 



^ (2z) fc fc!r(i/- k + \ 



2z 



+ exp 



z + 



is 2 - 1 



2z 



4 ±i7r(i/+i 



(3.24) 



("+" for -tt/2 < arg(z) < 3tt/2, "-" for -3tt/2 < arg(z) < tt/2). Langguth and 
Inomata [67] argued that the inserting of the expansion (3.20) in the path integral 
(3.17) can be justified (adopting an argument due to Nelson [80]), but things are not 
such as easy. Our arguments are as follows: 
1) The Schrodinger equation 



h 



d 2 D-ld\ r2 l(l + D-2) t r / , 
2m \dr 2 2r dr 2mr 2 { ' 



#(r;£) = iftJU(r;£) (3.25) 



can be derived from the short time kernel of the path integral (3.17). 
2) For r', r" — > equation (3.24) has the right boundary conditions at the origin. We 
have for z — > 0: 

V2^e~ z I x (z) ~ -JL-z*** [1 + 0(z 2 )] 



and it follows: 



(D) I s. 



(2A)! 

r a+(D-i)/2 (r / ^ q 5 r /, > Q) 

r ///+(D-l)/2 (/ / ^ ^ r , > Q) 



(3.26) 
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for all j. If and e^ denote the (normalized) radial state functions and 

energy levels, respectively, of the .D-dimensional problem, we have 



K\ >{r ,r-T) = ^R> nl , {r)R y n >{r )e 



JD)/n /h 



(3.27) 



n=0 



and thus for any but s-states - including the factor r ( D 1 )' 2 in equation (3.27) - 
we have: 

R (P) {r) „ rl+(D - 1)/ 2 r -(D-l)/2 = r l (3 _ 28) 

which is the correct boundary condition for R^i- 
Both points are not true for the path integral (3.22), because 

1) we did not find any obvious method for deriving the full Schroding er equation 
(3.25), the angular dependence comes out wrong. 

2) In the "Euclidean region" T — > — i t(t > 0) the functional measure is vanishing like 



exp 



<f + W - \ + (' + W - \ 



2mrir' 



2mr 




r N ^(l + 



D-2 \2 _ 1 
2 ) 4 



-n-Y 

N ^ 2mr (j) r (j _ 1) 



(3.29) 



which is also vanishing for r', r" 



0, but not in the right manner. The points 



r' = r = are essential singularities, where as the correct vanishing is powerlike. 
3) Equation (3.22) seems very suggestive, but it is quit useless in explicit calculations. 
Even in the simplest case D = 3, / = 0, V = 0, i.e. 



K 



l 



r",r';T) = (r'r") — ~ lim 



N^oo V27T 



m \ 



N/2 



X 



/oo roc \ im N 

dr (1 y ■ ■ dr^.i) exp I — ^(r 0) - r (i _ 1} ) 



(3.30) 



cannot be calculated. The integrals turn out to be repeated integrals over errorfunc- 
tions which are not tractable. Also the method of Arthurs [4] fails. In this method 
one assumes that the integrations in the limit e — > are effectively from — oo to +oo. 
The path integral (3.30) then becomes 



Ki(r",r';T) = (r'r") 2 exp 
which is wrong. The "mirror" term K^ D \r" ', —r';T) is missing. 



2eh 



(r' - r"f 



Let us discuss some properties of the path integral (3.22) 

r(t")=r" 



k 



(D) 



(r",r';T)= J ^ D) [r 2 ]Vr(t)e W Uj^ 



r(t')=r' 



y r 2 - V(r) 



dt 



(3.31) 
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1) The whole /- and D-dependence is contained in ^\°\ as defined in (3.18): 

d D) [z {j) ] = v^F^e-W) I l+ E^(z (j) ). (3.32) 
Therefore we can conclude: 



(3.33) 



and all dimensional dependence of the path integral (3.31) can be deduced from 
the three dimensional case: 



k ( l D) (r'\r'-T) = k^ £ ^(r'\r\T). 



(3.34) 



From now on we will denote fc,(T) = A; Z (3) (T). 
2) The boundary conditions we have already discussed; we have for the radial wave 
functions u\ D ^ and R\ D ^ respectively: 



uf\r) -> r l+£ ^ , R\ D \r)^r l (r -> 0). 



(3.35) 



We therefore state: The radial path integral fej(T) can be written as a superposition 
of two one dimensional path integrals: 



h{r"y;T) = ^(r^r'jT) - (-l^V, -r';T) 
with 



(3.36) 



jV /»00 /» 



dr 



(N-l) 



N ( . TV 



m 

27n 



A %)-^(|r (j) l) 



. (3.37) 



4) From equation (3.36) for the case / = V = we have two free particle path integrals 
the solution being: 



^=°(rV;T) 



m 



2n\hT 

m 
2ni HT 



exp ( — nrni __ (r — r 1 



exp 



2Th 
im 



— exp 



im 
2Tft 



(r" + r ; ) 2 



( r »_ r ')2 



2Th 



1 - exp ( - — r r 



(3.38) 

the second term is the so called "mirror" term. It is not allowed to drop it, even for 
very small T — > 0, since this violates the boundary condition which demands that 
kQ =0 vanishes for r', r" — > 0. These two parts of the Bessel functions express nothing 
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else than the famous mirror principle. Actually, in the above Euclidean form k$ is 
nothing else but the heat kernel on the half line r > 0, which has the representation: 



K (s) (r // j r >. T) = R (i) (r // j r >. T) _ K (i) _ r ,. T) 



(3.39) 



where K^ 1 ' denotes the one-dimensional heat kernel in R. In the language of diffusion 
processes: r = is an absorbing boundary and the job is done by the functional 

(D) 

measure n\ . 

Let us discuss the radial path integral in the language of the Weyl-ordering and 
product-ordering rule [49]. We consider the Schrodinger equation in D-dimensions 
with a potential V(|x|) = V(r): 



i ft J^*(r, {*};*) 



2m' 



+ V(r) 



(3.40) 



with Alb the Laplace Beltrami operator. Introducing the D-dimensional Legendre 
operator: 

d 2 



T 2 



d 2 d 
—2 + (D-2)cot6»i — 



d9 2 



+ 



1 



sin 9\ . . . sin #d-3 



d9 1 

d 2 



+ 



1 



dei_ 2 



sin 2 9 X 

+ COt 0D-2 



_ + (fl _3)cot^ 



+ 



d 



89 



D-2 



+ 



1 



d 2 



sin 2 6»i . . .sin 2 6>,d_2 <9(/> 2 



(3.41) 



the Laplace operator is: 



A LB = TT^ + 



8 2 D - 1 d 1 



dr 2 r dr r 2 ^ 



(3.42) 



Rewriting the Hamiltonian if = — t^Alb yields: 
H(p r , r, $}) = || 



+ 



1 



2mr 2 



2 1 2 1 



sin 2 X J 



sin 9i . . . sin 6 D -2 



with 



AV Weyl (r,{e}) 



h 2 



8mr 2 

and the hermitean momenta 



1 + 



1 



sin 2 0i 



+ ... + 



+ AV Weyl ({9}) (3.43) 



(3.44) 



sin d\ . . . sin Ob -2 







7 




\ 








ft 




i 90' 



+ 



cot 9 U 



(3.45) 
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Because of the special nature of g a b, the product ordering gives the same quantum 
potential as the Weyl-ordering, i.e. we have AVw ey i = AV pro d- In particular we have 
the equivalence 



N-l 



N 



II V^-exp ^^^^),^)},^- 1 ),^- 1 )}) 



N 



-{g'g")- 1 II V^-exp^^lf^,^')}) \ (3.46) 



where 



(flC?) _ ^- 1 ))2 + sin 2 ^^ - fl^) 2 + 



2/,0)//,0) _ /j(i-l)\2 



+ (shT 2 ^ j) . . .sm 2 ^:^)^ - ^- 1 ))^ 



(3.47) 



denotes a "classical Lagrangian" on the lattice. C^^f^', {9^}) is similar defined as 
(3.47), except that one has to take all trigonometries at mid-points. The features are, 
of course, direct consequences of the product form definition as described in section 
II. 3. With the correct Hamiltonian (3.43) we can consider the Hamiltonian path 
integral for the Feynman kernel K: 



K( D \r",{e"},r',{e'}-T) 

r(t")=r" 

{6{t")=e"} 

= WV* I Vr(t)Vp r (t)V{8}(t)V{p e }(t) 



x exp 



r(t')= 


1 

T 




0'} 


+" 




/ 


Pr 


Jt' 





D-l 



dt 



xe*p<- 



N 



J2Pr U) Ar U) + W)}A{^0)} - eH ef f(pr U) ,r {j) , {pe U) }, {%)}) 

(3.48) 
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with the effective Hamiltonian 



n effiPr u yr {j) ,{pe {jV e {j) }) 



p; 



2m 



+ 



2mrf j) 



i 



sin d\ J ' 



sin 2 ^...sin 2 ^l 2 9U) 



+ AV Prod ({6 {j) }) 
(3.49) 



After the integration over all momenta we get (AV = AVweyl = AV pro d): 



K^(r",r',{9"},{9'};T) 

^ 

Vr{t) DQ(t) exp 



r(t")=r" 
{8(t")=9"} 



h L, 



[C^(r,f,{9,9})-eAV(r,{9})}dt\, 



r(t')=r' 
{9{t>)=9>} 



m 



lim ( 

N^oo \2n\en 



ND N—l „oo 
2 



n 



'0') 



x exp 



i 



AT 



(3.50) 



with the effective Lagrangian 



£ eff (r,f,{9,9}) 



m , 



[r 2 + r 2 # 2 + r 2 sin 2 #^2 + • • • + r 2 (sin 2 9 1 . . . sin 2 D _ 2 )<^] - K(r, {6>}) 

(3.51) 

defined in the same way as in equation (3.47). 

The path integral (3.50) with the Lagrangian given by equation (3.47) is too 
complicated for explicit calculations. We therefore try to replace equation (3.47) 
under the path integral (3.50) by the following expression: 

~ • 777 

Cci(r, {#}, r, {9}) - C C i(r, {9}, r, {9}) := | R 2 Q 2 - V c ({9}) (3-52) 

where V c has to be determined and O denotes the D-dimensional unit vector on the 
S 0-1 -sphere. Thus we try to replace Cci by a simpler expression and hope that 
V c + AVweyl is simple enough. We have 



(0(D_0(2) ) 2 = 2(1 _ cosV ,(l,2) ) 
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with the addition theorem (3.3). We shall use equation (3.53) to justify the replace- 
ment (3.52) and thereby derive an expression for V c . We start with the kinetic term 
(x^ — x^ -1 ') 2 expressed in the polar coordinates (3.2), R = r (not fixed), and expand 
it in terms of Ar and A6 U . After some tedious calculations we obtain the following 
identity 



exp 



^C N cl {^\{e^r^\{6^}) 



exp 



lm 

eh 



( r (j)2 + r (j-i)2 _ 2r 0') r 0-i) cosTpti'i- 1 )) -i e v c (r^\ {6 ij) }) 



(3.54) 



with 



V c {r^\{9^}) 



8mrW 2 



1 + 



sin 2 ^ 



H + 



sin 2 #i j) ...sin 2 ^l 2 



(3.55) 



(V c is the same whether or not Ar^ = 0). The result is that the potential V c generated 
by these steps cancels exactly AVweyl(r, {0})\ Therefore we get: 



t")=r" 

,//\ n//~i 



{d(t")=e"} 



K^ D \r",r',{e"},{0 / };T) = J Vr(t) TXl(t) exp U J 



m 



-x 2 - V(r) dt 



r(t')=r' 

W)=o} 

(3.56) 

where x 2 has to be expressed in polar coordinates. In the lattice formulation x 2 reads 
• 2 ' [r 2 j} + r 2 ^) - 2r (j) r _ 1) cos iP UJ _ 1) ]/e 2 . (3.57) 



x 



Therefore we arrive at equation (3.52) and both approaches are equivalent as it should 
be. 

Let us note that equation (3.54) as derived in [49] is effectively a regularization 
of the highly singular terms in AV. This corresponds to the l/r 2 -terms in the radial 
path integral which have been regularized by the functional measure m [r 2 ] . However, 
the functional measure method gives a quite practical tool to regularize such singular 
terms. 



3.2. The Radial Harmonic Oscillator 

Let us now discuss the most important application of equation (3.17), namely the 
harmonic oscillator with V(r) = ^muj 2 r 2 . The calculation has first been performed 
by Peak and Inomata [84]. However, we present the more general case with time- 
dependent coefficients following Goovaerts [39]. This example will be of great virtue 
in the solution of various path integral problems. 
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We have to study 
tf,(r",r';r) 

= (r'r^l™ o (^ E ) N/2 l dr {1 y..J dr {N _ 1} 



N 



x n {A [ r o) r (j-i)] • ex p 

i=i L 

N/2 ,. x 



im o ie 

2^( r 0') _r (i-i)) 



2 ^m^(/'")r 



= (r r ) 2 lim , , 
tv^oo \ieh J 



uj I r(i)r(i) --y 



r (i) r (i)'" / r (J v_i)rfr (JV _i) 



N 



exp 



im 



ss (4, + ^- 1 ,)-i- 2 (' u, )4, 

/ \ N / 2 

lim 



•I,,n- 



m 



ir'r") 2 lim 



2 ia(r J )/2 



i / 



L 



r(i)dr {1 y ■ ■ / r (Ar _i ) <ir (A r_i ) 



x exp 



i{P(i)T 2 (X) + P(2)rf 2 ) + ■■■ + P(N- l)^(AT-l)) 



X 



^1+^=2 (- i « r (o)^(i)) • • • ^+£=2 (- i a^(7v-i)^(7V)) 



2 

2-D 



= (r'r")^" lim <(T) 



(3.58) 



where (T) is defined by the iterated integrals. Furthermore we have set a = raj eh 
and = a[l - e 2 muj 2 (t^)/2]. We are now using [40, p. 718] 

pe"^ 2 J„{ax)J„{l3x)dx = ±- e -(<* 2 +P a )/*i h (3 . 59) 

which is valid for > — 1, | arg ^fy] < ir/4 and a, (3 > 0. With analytic continuation 
[84] one can show that 



Jo 



re iar 7 1/ (-iar)/ I/ (-i6r)dr= -W a +& )/4ai ^ 

2a 



a6 
2ai 



(3.60) 



is valid for v > — 1 and K(a) > 0. By means of equation (3.60) we obtain for Kf 1 {T): 



K?(T) 



a 



exp I — (r + r 



POO POO 

/ r {1) dr {1 y / r (A r_i ) rfr (Ar _ 1) 
Jo Jo 



x exp 



' ; (Al) r ?i) + A2) r (2) + ' • ' + AiV- l)^(iV_l)) 

h+Zji ( _iar (o) r (i)) • • •^ + ^ 2 (- i «^(7v-i)^(JV)) 
— exp(ipArr +iQArr )I l+ n~2{-ia N r r ) 

1 2 
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X 
./2 
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where the coefficients cxn,Pn and qN are given by 

N-l 



fc=l 



a 



N-l 



Pn 



a 



a. 



1 ^ 4 v 

k=i /fc 



a a 
qN= 2~4^ 

ol\ = a, 



k 



a k +i = a J] 



a 
Wk 



(k>l) 



7i= A, 7fc+i=/? fc+ i 



(3.62) 



We must now determine these quantities. Let us start with the evaluation 
Putting 

27fe 2/fc+i 



we obtain 



7fc+i = Pk+i 



a 



47fc 



Vk+i - 2y k + y k -i 2 
o + u kVk = 0. 



of 7fc- 
(3.63) 

(3.64) 



In the limit N — > oo this gives a differential equation for y 

y + u; 2 (% = 

with solution y = + t) = 77 (t). Since on one hand side 

2/i - 2/o + e?/o -> 2/0 (e -> 0), 
2 

2/i - -2/o7i -»■ 2 2/o, (e->0), 
a 

and on the other 

2/1-2/0 1/2 \ 

2/0 = = - -7i-l 2/0^1, (e^l), 

e e \a / 

we have the boundary conditions 

77(0) = 0, 77(0) = 1. 
Now let us expand 2/fc+i — + l)e] + 0(e 3 ), we observe that 

Vk+i ~ 22/fc + Vk-i + e 2 LV 2 (t' + ke)y k = 
52 



(3.65) 
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is satisfied up to second order in e. Therefore expanding (3.64) 

= ayk+2 = a y[(k + 2)e}+0(e 3 ) 
7fc+1 2y k+1 2# + l)e]+0( f 3)- 

Consequently 

m % I r 1 rj(je) + 0(e 3 ) 



lim cx 7v = lim —If 



JV^oo eh + l) e ] + 0(e 3 ) 

m 77(0) 



= lim 



m 



JV^oo eft r/(eiV) ^(T) ' 



(3.66) 



Similarly 



m 



lim q N = lim 

tv^oo at^cx) Zen 



1 



ry[(iV-l)e]+0(6 3 ) 



= lim 



n^+oo 2h er][Ne] 
Finally we must calculates pn- We have 

(a N ~ X a 2 \ 
lim pn = lim ( — — > — — I 

iV^oo N^oo \ 2 4>J k / 

m /l Ar_1 
= — lim ( 

2% n^oo V e 



TjfeiV] + 0(e 3 ) 
m ?7[iVe] - 7/[f + (N — l)e] _ m?)(T) 



(3.67) 



— lim ( - 

m 



- r? 2 [(fc + l)e] 



eft 



ry2(t) 



lim 

2^(T) jv^oo 



r/(T) 



- V(T) 



dt 



m 



2hr](T) 



an 



(3.68) 



Furthermore we find 



lim £(e) = lim ^ - 77(e) 



lim 



77(0) + e7)(0) 



n^oo e 

77(0) = 1 
77(e) 



77 2 (t) 



lim £ (e) = lim 

e^O V ' e^O \ e 



77(e) 



7)(0) [77(0) + e7)(0)]-e 
" ~o e[77(0) + ci7(0)] 

and therefore £(rj) is satisfying the boundary conditions 

£(0) = 1, £(0) = 0. 
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£(£) is found to satisfy the differential equation 

e + u; 2 (t)£ = 0. 

Therefore we have 



P(T) = lim p N = 



m£(T) 



(3.70) 



(3.71) 



jv^oo 2hri(T) 

and we have finally for the path integral of the radial harmonic oscillator with time- 
dependent frequency: 



= (r'r") — 



2-d m 



exp 



In particular for u>(t) = ui = const.: 



2% \n(T) + r}{T) 



mr r 



i](t) = — sino;(t — t'), fj(t) = cosu(t — t') 



£(t) = cos u(t-t') 

which yields the radial path integral solution for the radial harmonic oscillator with 
time-independent frequency 



= (r'r") — . 



2-d muj 



iTismuT 



exp 



™, ,2 . ,2. „ 

(r + r ) cot ujT 



muJr ' r ' 

I, i D-2 . 

t+ 2 Vi ft sin u/T 



. (3.73) 



It is possible to get the same result, if one starts with the .D-dimensional path integral 
in Cartesian coordinates: 



x(t")=x" 



K(x", x'- r) = J Vx(t) exp J ^ J [x 2 - u 2 x 2 } dt 

x(t')=x' ^ 1 

insert for every dimension the solution of the one-dimensional oscillator 



(3.74) 



K{xlx' k ;T) 



2-KihsmuT 



exp 



i muj 
2h 



2 , /2 N , m 2x' k xl 



[x' k +x' k )cotuT - - 



sinu/T 



(3.75) 



and uses equations (3.8) and (3.11). 

The next step is to calculate with the help of equation (3.73) the energy-levels and 
state- functions. For this purpose we use the Hille-Hardy- formula [40, p. 1038] 



t -a/2 
1-t 



exp 



1 1 + t 

(x + y) 

2 V y 'l-t 



^ ' n=0 



t n nle -±(x+y) 

r(n + a + 1) 



(3.76) 
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With the substitution t = e 2lujI , x = mojr' z /ft and y = moor' 2 /ft in equation (3.73) 
we get finally: 

oo 

Ki{r'\r';T) = ^ e" {TEn/H R l N {r')R l N {r") (3.77) 

7V=0 

E N = uh[ N + ^ j (3.78) 



2mw r(^ + l) /"mw^ 1 ^ 



^ r D-2 r (iv+M) ^ ft 



exp 



ft 



ft 



(3.79) 



The path integral for the harmonic oscillator suggests a generalization in the index 
/. This will be very important in further applications. For this purpose we consider 
equation (3.17) with the nontrivial functional measure n\ D \r 2 ) = Hi[r 2 ]: 

r(t")=r" 



1 



Kfr", r'; T) = — / Vr(t)^[r 2 ] exp { - 



r(t')=r' 



ft 7 t / 



f r 2 - K(r) 



dt ^ (3.80) 



The functional measure corresponds to a potential V\ = ft 2 ^rj? in the Schrodinger 
equation. Assuming that we can analytically continue in / — ► A with 9ft(A) > —1, then 

we get for an arbitrary potential V\(r) = ft 2mr f naively inserted into the radial path 
integral 



r(t")=r" 



— r 2 - ft 2 
2 



\2 1 

A - 3 m 2 2 

^ r 

2mr 2 2 



(3.81) 



r(t')=T" 



This path integral must now be interpreted in terms of the functional measure in 
equation (3.18). Define (z = mr^.i)^)/ ieft): 



N 



(3.82) 



Then equation (3.81) must be interpreted with the help of equation (3.82) as 

r(t")=r" 



/ ^HU 



r(t')=r' 



m -2 + 2 

—r 2 - ft 2 



\2 1 

A - 3 m 2 2 

— ^ — —a; r 

2mr 2 2 



dt 



r(t")=r" 



r(t')=r' 

V r'r"mu> 
iftsinwT 



y Pr(t)^ A [r 2 ]exp ^ jf (r 2 -w 2 r 2 )dt 



exp 



ima; ( ,2 

~2h 



(r' + r' ) cot cuT 



( a 7 



mur'r" 



i ft sin cuT 



(3.83) 
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This important equation has been derived by Peak and Inomata [84] with the 
help of the three-dimensional radial harmonic oscillator, and by Duru [26] with the 
corresponding two-dimensional case. Duru considered the radial potential problem 



V(r) = -~ + 6r 5 



(3.84) 



\muj 2 it is 



with some numbers a and b. Identifying a = (h /2m) (A 2 — \) and b 
simple calculation to express equation (3.83) in terms of a and 6, which is omitted 
here. However, these authors did not discuss this path integral identity in terms of the 
functional measure language. That such a procedure is actually legitimate is beyond 
the scope of these notes. It was justified by Fischer, Leschke and Miiller [35] for the 
radial path integral and for the Poschl- Teller path integral as well (see below), where 
one also has to be careful with the appropriate Besselian functional measure. In the 
functional measure interpretation equation (3.83) can be found in references [49] and 
[92]. Equation (3.83) is very important in numerous applications. 

Let us note the free particle case. In the limit w->0we obtain in equation (3.73) 



K \ u=Q \r",r';T) 



(r'r") 



2-d m 



exp 



im , 2 ,2x 
(r + r ) 



2hT 



mr'r" 



ihT 



(r r ) 



2 J dp p exp ^ - ) J l+£ -2 (pr') J l+£ -2 (pr") 



(3.85) 



with wave-functions and energy spectrum 



2D % 2 2 

*p(r) = r^^p- J l+ D-2(pr), E p = —. 
The one-dimensional case gives (i.e. motion on the half-line) 



(3.86) 



[Y n 2 p 2 

* p (r) = r^/p Ji (pr) = \ — sinpr, E p = — — 
2 V Tvr 2m 



(3.87) 



However, there is an ambiguity in the boundary condition for r = for D = 1. The 
present case here, i.e. ^ p (0) = 0, corresponds to a specific self adjoint extension of the 
Hamiltonian H = —h 2 d 2 /dx 2 for functions \1/ G L 2 ( [0, oo)) on the half- line. 

Finally we calculate the energy dependent Green function for the radial harmonic 
oscillator. Making use of the integral representation of the previous section we obtain 



Gi(r",r'; E) 

e iET ' n K l (r",r';T)dt 

-W 



2Kuj ' 2 ' 



muj 2 



(3.88) 



56 



III.4 



Other Elementary Path Integrals 



From this representation we can recover by means of the expansion for the T-function 
the wave- functions of equation (3.79). 

The corresponding Green function for the free particle has the form 

2m 



%{r'r") 



D-2 
2 



1 7rm 

2 



1 / mE 



i V 2h z 



(r' + r" - \r" -r'\) 



1 mE 



i V 2h z 



(r' + r" + \ r "-r'\) 



(3.89) 



h(r'r") 



x J 



i+- 



--%(r' + r" -\r" -r'\) 

2h 2 



H 



(i) 



mE 



(r' + r" + \r" -r'\) 



(3.90) 



where use has been made of the integral representation 



dx 



x 



a + b\ r ( a — b 



exp I — px 



2x 



2x 



= 2I V 



(3.91) 



4. Other Elementary Path Integrals 



There are two further path integral solutions based on the SU(2) [5, 25, 60] and 
SU(1, 1) [5, 61] group path integration, respectively. The first yields the path integral 
identity for the solution of the Poschl- Teller potential according to 



K^ PT \x",x';T) 

x(t")=x" 



- i 



Vx(t) exp <| — 



x(t')=x' 



m 



-x 



h 2 fa 2 



+ 



P 2 



2m \ sin 2 x cos 2 x 



1=0 



exp 



ihT 

~2m 



(a + (3 + 21 + if 



(x')^ p> (x") 



= V sin 2x' sin 2x" exp 



j=o, 



ihT 
~2~m 



(2J + lf 



dt 



x (2J+ l)DiU i -^(cos2x')Di ±£ ^ 

2 ' 2 2 ' 2 



(cos2a;") (4.1) 



with the wave-functions given by 

V^(x) = y / (2J + l)sin2x Di 



<A j , p- a (cos 2x) 

2 ' 2 



(4.2a) 
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2{a + (3 + 21 + I) 



l\T(a + (3 + 1 + 1) 



x 



T(a + l + l)T(/3 + l + l) 
(smx) a+ ?(cosxf + ?Pi a > (3) (cos2x) 



(4.2b) 



a + (3 + 21 + 1 l\T(a + P + l + 1) 



2 a+/3+i r(a; + z + i)r(/? + Z + 1)_ 

x v / 2sin2a;(l -cos2a;)^(l + cos2x)^P^ a,/3) (cos2x). 

(4.2c) 

Here, of course we can analytically continue from integer values of m and n to, say, 
real numbers a and /?, respectively . 

Similarly we can state a path integral identity for the modified Posch- Teller 
potential which is defined as 



h fv 



2m V sinh r cosh r 



(4.3) 



This can be achieved by means of the path integration of the SU(1, 1) group manifold. 
One gets 



N M 



K( MPT )(rV';T) = ^$^ )t (r')$^ ) (r")exp - -f- 2(h - k 2 - n) - 1 

**' ill 



n=0 



ihT 



poo 

+ / dp^' u) *(r')^' u) (r")exp 
Jo 



ihT 
~2m 



V 



(4.4) 



Introduce the numbers fci, k 2 defined by: k± = |(1 ± u), k 2 = \{1 ± rj), where the 
correct sign depends on the boundary conditions for r — > and r — > oo, respectively. 
In particular for rf = |, i.e. /c 2 = \, §, we obtain wave-functions with even and odd 
parity, respectively. The number A/m denotes the maximal number of states with 
0, 1, . . . , N M < ki — k 2 — \. The bound state wave- functions read as: 



$(fci.fc2)( r ) = ^' fc2 )(sinhr) 2fc2 -5(coshr)- 2fcl+ i 

x 2Fi(-k 1 + k 2 + k, -ki + k 2 - k + 1; 2k 2 \ -sinh 2 r) 



jy(fci,fc 2 ) _ 



r(2A; 2 ) 



2{2k - l)r(/ci + fc 2 - Av)r(/ci + k 2 + k — 1) 
r(fci - fc 2 + /c)r(fci - fc 2 - « + 1) 



> (4-5) 



(k = k± — k 2 — n) , (there is a factor "2" missing in reference [36]). Note the equivalent 
formulation 



$( r ,/3) (r) = 



(3n\T(n + a + (3 + 1) 
2 Q +/ 3 r(n + a + l)T(n + 13 + 1) 



(l-x)^(l + x)^Pi a ^(x), (4.6) 
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with the substitutions a = 2k 2 — 1, (3 = 2(ki — k 2 — n) — 1 = 2k— 1, x = 2/ cosh 2 r— 1 with 

r -i -l 

the incorporation of the appropriate measure term, i.e. dr = (1 + x) a/2(1 — x) dx. 
The scattering states are given by: 



$(fci.*a)( r ) = ^ fcl ' fc2) (coshr) 2fcl -5(sinhr) 2fc2 -5 



1 



fj{ki,k 2 ) — 

p T(2k 2 ) V 2tt 2 



x 2-Pi(&i + fc 2 - «, fci + fc 2 + « - 1; 2/c 2 ; - sinh 2 r) 

r(fei + fc 2 - «)r(-fei + fc 2 + «) 



p sinh up 



(4.7) 



x r(fei + k 2 + k - i)r(-fei +k 2 -n + i) 



[K=\{l+ip)}. 

It is possible to state closed expressions for the (energy dependent) Green functions 
for the Poschl- Teller and modified Poschl- Teller potential, respectively. For the Poschl- 
Teller potential it has the form (Kleinert and Mustapic [65]) 

G(x»,x';E) = ^V sin2^sin2, % , r(mi = + ™* + *) 

v ; ih r(mi +m 2 + l)r(mi - m 2 + 1) 



x 



X 



1 - cos 2x' \ ^- m *V 2 ( 1 + cos 2x' \ ( mi + m2 )/ 2 



1 - cos 2a:" \ (mi " m a)/ 2 ( 1 + cos 2x" \ ^ +m ^/ 2 

1 — cos 2x' 



x 2 Fi - L s + mi; L s + mi + l;mi - m 2 + 1; 



x 2 F 1 - L E + mi; L E + mi + 1; mi + m 2 + 1; 



1 + cos 2x" 



(4.8) 



with mi/ 2 = | (A ± k), L s = — | + \\j2mE /% and x" > x'. A similar expression 
is valid for the modified Poschl- Teller-potential Green function (see [65]) which is, 
however, omitted here. 



5. The Coulomb Potential 

The hydrogen atom is, of course, one of the most interesting subjects in quantum 
mechanics. In the beginnings of quantum mechanics and with the atom model of 
Rutherford it was a riddle how to tract a system which is from the classical physics 
point of view unstable and doomed to vanish into pure radiation. Its was Bohr's 
genius, postulating the famous rules that only a countable number of orbits are allowed 
satisfying the quantum condition § pdq = nh (n G N). However, this "old" quantum 
mechanics was not sufficient, because it fails e.g. in the case of the helium-atom, and 
it was Pauli who solved the hydrogen problem in the terms of the "new" quantum 
mechanics developed by Schrodinger and Heisenberg. It is surprising, that Pauli did 
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not use a differential equation and solves the corresponding Eigenvalue problem as 
Schrodinger did later, instead he exploited in fact the "hidden" SO (4) symmetry of 
the Kepler- Coulomb-problem. This symmetry gives classically rise to an conserved 
quantity, called Lenz-Runge vector. This additional symmetry allows also a separation 
of the Coulomb problem in parabolic coordinates. This vector points from the focal 
point of the orbit the perihel. 

Ever since the success of quantum mechanics the hydrogen atom was the model 
to test the theory, may it be Dirac's relativistic quantum mechanics, where first the 
fine-structure constant a = e 2 /he arises. 

It was for a long time a really nuisance that this important physical system 
could not be treated by path integrals. Calculating wave functions and energy 
levels remains more or less a simple task in the operator language, but even the 
construction of the Green function (resolvent kernel) was impossible for a long time. 
It takes as long as 1979 as Duru and Kleinert [27, 28] finally applied a long-known 
transformation in astronomy (Kustaanheimo-Stiefel transformation [66]) in the path 
integral of the Coulomb problem and were successful. Even more, their idea of 
transforming simultaneously coordinates and the time-slicing opens new possibilities in 
solving the huge amount of unsolved path integral problems. In particular, Coulomb- 
related potentials, in the sense, that all these problems can be reformulated in terms 
of confluent hypergeoemtric functions, like the Morse-potential could be successfully 
treated. However, the original attempt of Duru and Kleinert was done in a more or 
less formal manner, and it does not takes a long time when it was refined by Inomata 
[58] and Duru and Kleinert [28]. 

Closely related to the original Coulomb problem is, of course, the 1 /r-potential 
discussion in R D . The D = 2 case in discussed in this subsection, whereas the original 
Coulomb problem in the next. 

As it turns out its Schrodinger equation for the Coulomb potential is separable in 
four coordinate systems: 

1) Spherical coordinates: 

x = r sin 9 cos <p 

y = rsm9sm<f) (r > 0, < 9 < tt, < </> < 2tt). (5.1) 
z = r cos 9 

Separation of variables in this coordinate system is, of course, not a specific feature 
of the Coulomb problem, but a common property of all radial potentials. 

2) Parabolic coordinates: 

x = £r] cos (j> 

y = £r7sin0 (£,??> 0,0 <</><27r). (5.2) 

Here electric and magnetic fields can be introduced without spoiling separability. 
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3) Spheroidal coordinates: 

x = sinh £ sin rj sin <p 

y = sinh £ sin 77 cos (£ > 0, < i] < tv, < <fi < 2%) . (5-3) 
z = cosh £ cos rj + 1 

Here a further charge can be introduced and therefore these coordinates are 
suitable for the study of the hydrogen ion H^- 

4) Spheroconical coordinates: 

x = rsn(a,k)dn(P,k') r > 0, \a\ < K, \/3\ < 2K', 
y = rcn(a,k)cn(p,k') 2 (5.4) 



z = r dn(a, k) sn(/3, k') 



0<k,k'<l, k 2 + k' =1. 



sn(ct, /c), cn(a, k) and dn(a, /c) denote Jacobi elliptic functions (e.g. [40, pp.910]) 
of modulus k and with real and imaginary periods 4K and 4iK' , respectively. 
Here the corresponding wave functions in a and (3 can be identified with the wave 
functions of a quantum mechanically asymmetric top. 

We will discuss the path integral for the Coulomb system in the first two of these 
coordinate systems. For the usual polar- and parabolic coordinate coordinate system 
the path integration can be exactly performed. In the remaining two, however, the 
theory of special functions of these coordinate is poorly developed and no solution 
seems up to now available. Nevertheless it is possible to formulate the Coulomb- 
problem in these coordinate systems and point out some relations to other problems 
connected to these coordinates [46]. 

It is surprising that the Coulomb path integral was first solved in Cartesian 
coordinates (where the Kustaanheimo-Stiefel transformation works) and not in polar 
coordinates. Looking at the appropriate formulas we see that even in the one- 
dimensional case (in polar coordinates) we need the one- dimensional realization of 
the Kustaanheimo-Stiefel transformation and a space-time transformation. 



5.1. The 1/r-Potential in R 2 [28, 57] 

We consider the Euclidean two-dimensional space with the singular potential V(r) = 
—Ze 2 /r (r = x G R 2 ). Here, e 2 denotes the square of an electric charged Z 
its multiplicity (including sign). However, as already noted, this potential is not the 
potential of a point charge in R 2 . The classical Lagrangian now has the form 



C(x,x) = —x + 



(5.5) 



Of course, bound and continuous states can exist, depending on the sign of Z. The 
path integral is given by 



x(t")=x" 



K(x",x';T)= J Vx(t)exp 1 jf ^ + ^f) 



dt 



x(t')=x' 



(5.6) 
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Discussion of this path integral are due to Duru and Kleinert [27,28] and Inomata 
[58]. However, the lattice-formulation is not trivial for the 1/r-term. In fact, it is too 
singular for a path integral, respectively, a stochastic process and some regularization 
must be found. This is known for some time, and it turns out that the Kustaanheimo- 
Stiefel transformation does the job. 

As it turns out one must perform a space-time transformation in order to solve 
the path integral (5.6). We perform first the transformation 

Xl = i 2 - r/ 2 , x 2 = 2£77, u = ( G R 2 , (5.7) 

which casts the original Lagrangian into the form 

C(u, u) = 2m(e + V 2 )(i 2 + V 2 ) + 7^-2 • ( 5 - 8 ) 
The metric tensor (g a b) an d it's inverse (g ab ) are given by 



1 0\ , ah , 1 (I 



with determinant g = det(g ab ) = 16(£ 2 + if) 2 and dV = dx x dx 2 = 4(£ 2 + if)d^drj. 
The hermitean momenta corresponding to the scalar product 

/OO /* OO 

/ d^ie+v^iitvm^v) (s.io) 
-oo J — oo 

are given by 

pc = tU + ^t^J' ^ = tU + ?t^J- (5 - n) 

Following our general theory of Chapter II. 5 we start by considering the Legendre 
transformed Hamiltonian 

h 2 ( d 2 d 2 \ ze 2 , s 

h ° = -^{m + m)-—- k (5 ' 12) 

which gives 

Therefore the time-transformation is given by e = f(£,rj)8 with f(£,r)) = 4(£ 2 + rf) 
and the space-time transformed Hamiltonian has the form 



*=-Uw + w)-* Ze2 -* E{eW) 

= 7^(P 2 + V%) - 4Ze 2 -AE(e + rf) 
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with the momentum operators and vanishing quantum potential AV: 

h d h d 



PC = 



id? 



i dr/ 



AV = 0. 



(5.15) 



This two-dimensional transformation can be interpreted as a two-dimensional Kustaan- 
heimo-Stiefel transformation [66]. It has the general feature that it performs a change 
of variables to "square-root" coordinates, note r = £ 2 + r] 2 . The general properties of 
transformations like this states the following problem (see e.g. [28, 46] for some review 
of the relevant literature) : For which values of D does exist a formula 



(xl + --- + x 2 D )(yl + --- + y 2 D ) 



4 + 



+ 4, 



(5.16) 



where the Zi are homogeneous bilinear forms in x and y. Following a theorem of 
Hurwitz and Lam this type of transformation can only be realized in the space- 
dimensions D = 1, 2, 4, 8. The assumption about z now implies that 



z = B{x) • y 

z i = J2(B k ) l i x k y l = J2B(x) l iyi- 

k,i i 



(5.17) 
(5.18) 



For the special case x = y one has 



x. 



(5.19) 



and the matrix B satisfies the condition 



t B(x)B(x) = \x\ 



(5.20) 



The one-dimensional case is "trivial", but we shall use it in the calculation of the 
path integral for the 1/r-potential in polar coordinates. The D — 2 case we have just 
encountered. The four-dimensional variation of this transformation has been used a 
long time ago in astronomy for the purpose of regularizing the Kepler problem. The 
D = 8 case is degenerate. Here the matrix B has the form: 



B{x) 



( 


X\ 


x 2 


x 3 


X4 


x 5 


Xq 


x 7 


x 8 \ 




-x 2 


X\ 


X4 


-x 3 


Xq 


-x 5 


x 8 


— x 7 




-x 3 


— X4 


X! 


x 2 


X 7 


-x 8 


-x 5 


Xq 




—X4 


X 3 


-X2 


X! 


-x 8 


-x 7 


Xq 


X 5 




-x 5 


-Xq 


-x 7 


x 8 


x 1 


x 2 


X 3 


— X4 




-Xq 


X 5 


x 8 


X-j 


-x 2 


X\ 


—X4 


-x 3 


\ 


—X7 


-x 8 


x 5 


-Xq 


—x 3 




Xi 


x 2 


-x 8 


x 7 


-Xq 


-x 5 


X4 


x 3 


—x 2 


Xi / 



(5.21) 



and maps R — > R and is thus of no further use. 



63 



Important Examples 



In the next Section (hydrogen-atom) we need the D = 4 case which is more involved 
and not one-to-one. 

To incorporate the time transformation 



8(t) 



1 da 



s" = s{t") 



(5.22) 



and its lattice definition e — > Af^'As^' = Au^' 2 8^' into the path integral (5.6) we 
now use from the set of equations (II. 4. 35), namely equation (II. 4. 35a). Observing 
that in the path integral (5.6) the measure changes according to 



N 

n 



m 



N-l N , . N-l 

x n = i n \^hm) - n ^ 



7 j=i i=i v 

we arrive at the path integral transformations equations (u G R 2 ): 

1 /"°° 

K(x",x';T) = — — / dE e- iTE / h G(x",x';E) 

2TTlh J_ QO 
nOO 

G(x'\ x'- E) = i / cfc" u'; s") + ^(-u", u'; s") 

Jo L 

where the space-time transformed path integral K is given by 



(5.24) 
(5.25) 



xexp^- 



ds } . (5.26) 



Note that the factor 1/r" has been exactly canceled by means of equation (II. 4. 35a). 
This is a specific feature of this potential problem. Furthermore we have taken into 
account that our mapping is of the "square-root" type which gives rise to a sign 
ambiguity. "Thus, if one considers all paths in the complex x = x\ + «X2-plane from 
x' to x", they will be mapped into two different classes of paths in the w-plane: Those 
which go from u' to u" and those going from u' to — u" . In the cut complex x-plane 
for the function u = Vx? these are the paths passing an even or odd number of times 
through the square root from x = and x = — oo. We may choose the u' corresponding 
to the initial x' to lie on the first sheet (i.e. in the right half w-plane). The final u" 
can be in the right as well as the left half-plane and all paths on the x-plane go over 
into paths from v! to u" and those from u' to — u" [28]". Thus the two contributions 
arise in equation (5.25). Equation (5.26) is interpreted as the path integral of a two- 
dimensional isotropic harmonic oscillator with frequency to = ^—8E/m. Therefore 
we get 



K{u", u'; s") 



2 i h sinews 



77 ex P 



4i Ze 2 s" 

h 



muj 
2nih 



(u' 2 + u' 2 )cotus" -2 



u' ■ u" 



sm us' 



(5.27) 
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Introducing now two-dimensional polar-coordinates 



£ = v /Fcos^, r7 = v^sin^, (r > 0, 6 e [0, 2tt]) 



so that u' • u" = y/r^r" cos(6" - 0)/2. Thus 
G(x",x';E) ' 



Tvfl J Q Sill ^'.S' 



x exp 



4iZe 2 .. moo , . .., 

— s -2Ti (r + '- ) 



, , mooy/r r" 6—6 

cosh ( — cos 

lfismws" 2 



Using the expansion 



cos 



(*«»£)= E 
we get for the kernel G(x", x';E): 

G(xf',xf;E) = -±- Gi(r",r';E), 



(5.28) 



(5.29) 



(5.30) 



(5.31) 



where the radial kernel Gi(E) is given by 
Gi(r",r';E) 



Imoo 



ds" 



sinws 



77 exp 



4iZe 2 s" moo , 



h 



-^-(r +r") cot ujs" 
a \ n 



<-2l 



mooyr'r" 



i h sinks' 



(Set oo = substitution u = and Wick-rotation) 



2m 



du 



h J sinhw 



exp 



2i 

ap 



u + \p{r' + r") cothw 



/ 2pV r'r" 
21 \ i sinhw 



(Substitution sinh-u = 1/ sinh-u) 

2 

2m f°° ( v\°^ 
= J ^coth-J exp [ i p(r' + r") cosh v\ In ( — 2 i p\J r'r" sinh v) dv 



(Reinsert hp = \JlmE) 
1 



m 1 /l , Ze 2 



VrV'V 2S(2/)! V 2 

ft V 2E' 1 



m 
2~E 



8m£ 



v- 



8mE 



(5.32) 



In the last step we have used (2.22) and r>,r< denotes the larger (smaller) of r',r", 
respectively. The representation (5.32) shows that Gi(E) has in the complex energy- 
plane poles which are given at the negative integers ni = 0, —1, —2, ... at the argument 
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of the r-function and a cut on the real axis with a branch point at E = 0. Thus we 
get a discrete and continuous spectrum, respectively: 



En 



mZ 2 e 4 



2% 2 {N - i) 2 ' 



(N = l,2,...) 



n 2 v 2 



(5.33) 
(5.34) 



Here we have introduced the principle quantum number N := n\ + I + 1. This is the 
well-known result. 

To determine the discrete spectrum we consider the first line in equation (5.32), 
identify in the Hille-Hardy formula t = e~ 2u , x = mur' and y = mur" and thus get 



n=0 2 ap 



X 



(2l + n)\ 



exp[-ip(r' + r")]L {2l \-2ipr')L {2l \-2ipr"). (5.35) 



Taking equation (5.35) and the n th residuum gives the energy levels and the wave 
functions. Inserting into equation (5.29) thus yields the Green functions for the 
discrete levels for the two-dimensional 1/r-potential 



G(x",x';E) = hJ2 



oo oo 



JV=1 l=-oo 



En — E 



(5.36) 



and the wave functions of the discrete spectrum are given by: 



*iv,i(r, 4>) 



(N-l-l)\ 



ira 2 (N - f) 3 (iV + Z-l)! 



x exp 



2r 



a(iV-i) 



a(N-±) 



-(20 
J N-l-l 



2r 



a(iV-i) 



(5.37) 



Furthermore we have used the abbreviation a = T% 2 / (mZ e 2 ) (the first "Bohr-radius"). 
The correct normalization to unity can be seen from the relation 



x 2l+\+2 e -x 



,(2l+A+l), x 
^N-l-l \ X > 



dx 



2(N+±)(N + l + \)\ 
(N-l-l)\ 



(5.38a) 



respectively 



2 dx = (2n + // + l)r(n + /x + l) 



(5.38b) 
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In order to determine the continuous wave functions we use the dispersion relation 
[40, p.987] 



/oo 
dx e~ 2lxp Y{\ + v + ix)Y(\ + v - ix)M ix , u {a)M ix , u {p) 
-oo 



, exp - (a + /?) tanh p J 2l/ -^JL . 5.39 
cosh p L V cosn P / 



and get (%> = \j2mE ) 
G,(r",r';£) 



2ip 



ny/r'r" Jo 



exp 
sinews" 

<is" exp 



|r(| + / + ip) | 2 
r 2 (2/ + i) 



4iZe s 
ft 

4iZe 2 
ft 



exp 



s" -p(r' + r") cot us' 



/ 2py/r'r" 

hi — - 

V sm us" 



4 i whs" \ 

^ +np)M iPjl (-2ipr')M_ iPtl (2ipr")dp 



m 



= ft 



dp 



( 



o ■ exp - 

hp 2 /2m-E \ap 

r(± + i+ ±)r(± + i- ±) 

V2 ap I V2 00/ 



X 



2nyfi J V"(2l)\ 



M_^ iZ (2ipr / )M^ iZ (-2ipr"). 

n.r> ' n.n ' 



(5.40) 



Here again the residuum at E = h 2 p 2 /2m has been taken in the integral. Thus the 
wave functions of the continuous spectrum are given by (hp = V2mE): 



* P ,i(r, </>) 
Note that 

r 



i(27)! r G +/ + i )eXP 



7T 

2ap 



-il<f> M_i_,(-2ipr). (5.41) 



1 7 i 
- + /+ — 
2 



Mj_ ,(-2ipr) 

17. T) ' 



2 op 



M_^,(2ipr) (5.42) 



and the wave- functions are therefore basically real. These results of the 1/r potential 
for the discrete and continuous spectrum, respectively, are equivalent with the results 
of the operator approach. The complete Feynman kernel therefore reads 



00 00 



K(x", x>- T) = J2H *- iTEN/h (r ; , 00* N,i(r", <t>") 

1=0 N=l 

OO /»CXD 

+ £/ *e- IT ^/ R *; ii (r',0 / )*^(r /, ,0") (5-43) 

with wave functions as given in equations (5.37) and (5.41) and energy spectrum 
equations (5.33, 5.34). 
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5.2. The 1/r-Potential in R 3 - The Hydrogen Atom [27, 28, 41, 55, 57, 58] 

We consider the Euclidean three-dimensional space with the singular potential V(r) = 
—Ze 2 /r (r = \x\, x G R 3 ) with Z and e 2 as in the previous Section. This potential 
corresponds in the space R 3 , of course, to the potential of a point charge and is thus 
the kernel of the Poisson equation in R 3 . The classical Lagrangian has the form 



m 



C(x, x) = —x* + 



Ze 2 



The path integral is given by 

x(t")=x" 



K(x",x';T) = j Vx(t)exp (^x 2 + ^pjdt 



x(t')=x' 



(5.44) 



(5.45) 



This path integral was first successfully solved (actually the Green function or resolvent 
kernel) by Duru and Kleinert [27], followed by further contributions of Inomata [58] 
and Duru and Kleinert [28], Ho and Inomata [55] Grinberg, Maranon and Vucetich 
[41,42] and Kleinert [62]. 

The transformation in the two-dimensional case corresponds to the two-dimensional 
realization of the Kustaanheimo-Stiefel transformation and we must look for a general- 
ization. However, we have a the relevant space R 3 whereas the relevant Kustaanheimo- 
Stiefel transformation maps R 4 -> R 3 which means that the transformation in ques- 
tion is not one-to-one. We have namely [28] 




I u 3 u 4 

-u 2 -ui 

-U\ U2 

u 4 -u 3 



Ul 
U 2 



Us 
-U4 



\ 


/«1\ 




u 2 




us 


) 


\ U4 / 



(5.46) 



or in matrix notation x = A ■ u with x G R 3 and u G R 4 , respectively. Several 
modifications of the matrix A are used in the literature, e.g. [10,11] 



A = 



( Us 
—U4 
-ui 
-u 2 



U4 

Us 
-1*2 
Ul 



Ui 

u 2 

Us 

—U4 



u 2 \ 

-Ul 

U4 
-Us I 



(5.46b) 



and see [66] as well. The transformation has the property A 1 A = u\ + u\ + u\ + u\ 
r = \x\ and no Jacobean exist. However, a fourth coordinate can be described by 



£4(5) = 2 / (u4u[ — usu' 2 + u 2 u' 3 — uiu' 4 )da. 

Js> 



(5.47) 



Therefore we must circumvent this problem. The idea [27,58] goes at follows: We 
consider the lattice formulation of equation (5.45) and introduce a factor "one" by 
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means of 
1 = 



m 



2n\hT 



d£" exp 



m 



2\KT 



= lim 



m 



N 



n^qo \2Tvieh 
which gives for equation (5.45) 

K(x",x';T) = 



HI d^exp 



lm 

We 



(5.48) 



X 



i ^ 

exp UE 



i=i 



g(AV« + A^>) + £ §f 



(5.49) 



with the to-be-defined quantity r^^. We now realize the transformation A (it) on 

4 



midpoints [58] iti J ' = \{ua +Ua ) (we use in the following the abbreviation x 4 = £): 



Ax^ = 2^A ab (^' ) )A^' ) 



6=1 



We have 



A\u( j) )-A(u( j) ) = u¥ )2 + uZ>" + u%>" + u%^ =: r 



which defines Furthermore 



(5.50) 



(5.51) 



A 2 ^ + A 2 4 J) + A 2 ^ + A 2 xf 



= 4f ^ [A 2 ^ j) + A 2 u ( 2 j) + A 2 4 J) + A 2 4 J) ] =: 4f V>. (5.52) 

Infinitesimal this has the form dx^ = 2A ab (u^)dua (sums over repeated indices 
understood) and the Jacobean of this transformation exists and is given by 



d(x 1 ,X 2 ,X 3 ,£) _ 2 4 f (j)2 



d(u 1 ,u 2 ,u 3 ,u 4: ) 
Thus we arrive at the transformed classical Lagrangian 

Ze 2 



(5.53) 



C(u, u) = 2mu u + 



u 



2 " 



(5.54) 



We now repeat the reasoning of the previous Section by observing that under the 
time-transformation equation (5.22) 



f* da 



s" = s(t") 
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and its lattice definition e — > 4r^As^ = 4u^ 2 5^ the measure in the path integral 
transforms according to 



N 



N-l 



N 



Respecting again equation (II. 4. 35a) which produces a factor 4r" we get the path 
integral transformation 



K(x",x';T) = y dE e- iTE / h G(x" 7 x';E) 

/•OO 

G(x",x';E) = i / d S "e 4iZe2a "/ R Jr( u ", u ' ;s ") 



(5.56) 



where the space-time transformed path integral K is given by 



K(u", u'- s") 

2 r°° 



d£" lim 



4r" ./_oo "* oo\27riM 



27V 



AT-1 „ 

2ivih sin uis" 



exp 



• N / 

i=i v 



00 <£" 



exp 



(w/ 2 + u' 2 ) cot us" — 2 — 



it' • it" 



sinws' 



(5.57) 



Here we have used once again the known solution for the harmonic oscillator, where 
we have actually a four-dimensional isotropic harmonic oscillator with frequency 
uj = ^/—8E/m. Note that the factor 1/r" has not been canceled as in the two- 
dimensional case. 

Up to now we have not discussed the problem of an eventually quantum correction 
appearing in the transformation procedure. As usual we consider the Legendre 
transformed Hamiltonian 



h 2 Ze 2 ^ 
H E = -- — A ( 3 ) E. 

2m y ' r 



(5.58) 



Let us write He in the coordinates: 



u\ = u cos a cos (3 
U2 = u cos a sin j3 
U3 = u sin a cos 7 
U4 = u sin a sin 7 



(u = \u\ = y/r, a= |, (3 ± 7 = 0) 
(0 < 9 < 7T, < 4> < 2tt) 



(5.59) 
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Then the Schrodinger equation in polar coordinates 



h 2 f d 2 2d 1 r2 \ Ze 2 ^ 
^ + -^ + ^L 2 ) E 



2 \ dr 2 r dr 



*(r,(9,0) = O 



(5.60) 



is transformed into 



h 2 1 ( d 2 3d ~\ Ze 2 E1 " 
+ -— -AK 2 — -E 



with 



2 4s 2 \ds 2 s ds 
d 2 



#(u) = 



d 



da 2 



(tana — cot ot) — 



1 8 2 



cos 2 a 



d 2 



(5.61) 



(5.62) 



da sin 2 a dp 2 d^ 2 

is the Casimir-operator of the group SO (4). Writing back into iti, . . . , 114 we get for 
He in the coordinates q G R 4 : 



H = Au 2 H E = 



2m 



A 4 -4Ze 2 -4Eu z 



and with p Uk = — i hd/dui-'- 



! 4 



(5.63) 



fc=i 



and no quantum correction appears. Therefore equation (5.57) is correct. 

To perform the ^"-integration in equation (5.57) we introduce polar coordinates 



/- 6 a + 4> 
u\ = v r sin - cos 

Zi 



r ■ u ■ 
U2 = yr sin - sin 

Z 



u 3 = \A* COS - COS 

Zi 



= y/r cos - sin 

Z z ■ 



2 




a + 




2 




a — 


<P 


2 




a — 






/0<9<n \ 

0<(p<2n 
\0<a<AnJ 



(5.64) 



We have with equation (5.47) 

f(s) = - / (a'-cos00>(<7)d<7. 
./o 

The polar coordinates give the expansion 



(5.65) 



, 0' . 0" 



X 



sin — sin — cos 
2 2 



a" - a' + 6" - 



9' 6" 
+ cos — cos — COS 
2 2 



a' - a" + 6" - 



(5.66) 
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Using the expansion 



e zcos *= e 1 (5-67) 

v= — oo 

twice, respecting due to the explicit form of £(s) from equation (5.65) above 

r°° d^) = r~ 

J-oo r " JO 



da 



(5.68) 



and that the integration over a" yields 4ir8„ ll/2 we get for the resolvent kernel: 



G(x", x>- E)= l - ds" - exp 

tv J \ 2 i n sin us" J 

OO / 



AiZe 2 s" mu 



h 



-^-(r +r") cot us" 
Zt \ n 



mu^fT? 1 . 0' . 9" 



sin — sin — — — 
hsinus" 2 2 1 ' p ' 



muVr'r" 9' 9" 



cos — cos — 
in sinus 2 2 



(5.69) 



The last sum can be exactly performed with the help of 



zi - z 2 e 



z\ - z 2 e 1 ' 
in particular 



I v Uz 2 1 +Z 2 2 -2z 1 Z 2 COS<Pj = J2 (-l) n In(z2)Iu+n(zi)e in( f > , 

^ ' n= — oo 

(5.70) 

j>2 j"* I v {z)I v {z') = I (\/z 2 + z' 2 + 2zz'cosA , (5.71) 

v=— oo ^ ' 

and the trigonometric expressions sin 2 ^ = ^(1 — cosa) and cos 2 ^ = |(1 + cosct). 
Thus 



G(x", x'- E) = - ds" 

TV 



mu 



2 i h sinus" 



x exp 



4 i Z e 2 s" mu , , 



h 



^ir' + r") cot us" 

2ih 



muy/r r ' 7 

cos — 

in sinus 2 



(5.72) 



where cos 7 = cos 9' cos 9" + sin^'sin^" cos(<p" — (p'). The Jo-Bessel function can be 
expanded in terms of Legendre-polynomials: 



(\ 2 °° 
zcos ^) =-^(2l + l)P|(cos 7 )J 2i+1 ( 2 ). 
l ? z 1=0 



(5.73) 



This relation can be derived from the general expansion 

7 \ /I V OO TV I 7\ 

(5.74) 
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Equation (5.73) is recovered with the identifications /j, = 1, v = 0, k = cos ^ and 
Pi(z) = (—l) l 2Fi(—lJ + l;l;^p). Thus the resolvent kernel can be expanded 
according to 



G(x", x'; E) = G(r", 9", <j>", r', 0', (/>'• E) 

oo 

= — J>z + i)PKcos 7 )G/(r"y;£) 



47T 

1=0 

oo I 



= E E Y l n *(e',<P')Yr(e",<l ) ")G l (r",r';E), 



1=0 n=-l 



(5.75) 



(5.76) 



where the radial Green function Gi (E) is given by 



G^'y-E) 

2mw 



00 ds" 



HVr'r" Jo sincus 
1 



exp 



4 i Z e 2 s" mw 



ft 



2ift 



(r' + r") cot us" 



<2Z+1 



y/r'r" 



iftsincus" 



rV'V 2£(2Z + 1)! 



m 1 „ / , i Z e 2 / m 



h V 2E 



, -8mE , , , 

X W Zc 2 /—xr,,i \ — -o — r> M Zo 2 j 

^fi - V _ 2B' <+ 2 \ V ft I h V ~2E ' ' 2~ 



=^r<|. (5.77) 



From the poles of the T-function at z = —n r = 0, —1, —2, ... we read off the bound 
state energy levels which are given by 



En 



mZ 2 e 4 
2ft 2 iV 2 ' 



( N = n r + 1 + 1 = 1,2,3,...) 



(5.78) 



which are the well-known Balmer levels. 

The determination of the wave functions is analogous to the previous two- 
dimensional case and thus we only state the result for the bound state wave 
functions of the hydrogen atom: 



(N-i-iy. 

a 3 {N + l)\ 



exp 



(5.79) 

A quite complicated attempt was made by Grinberg, Maranon and Vucetich [41] 
to determine these functions from the wave functions of the corresponding four- 
dimensional oscillator. The continuous state wave functions of the hydrogen 
atom read as: 



- r( 1 + / + — ] exp ( — )Mi i(-2ipr)K n (6>,0). 

27r(2/ + l)!r V <W \ 2 <*P J ^' +5 Ml ' 

(5.80) 
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Of course, these wave functions form a complete set. The spectrum of the continuous 
states is, of course, given by 

E p = (p > 0). (5.81) 

These results coincides with the operator approach. The complete Feynman kernel 
therefore reads 

oo oo 

K(x", x'; T)=J2J2 e~ iTEN/h ** N ,i(r', 0')*Ar,,(r", 0") 

1=0 N=l 

+ E/ dpe- {TE *l n %j{r'^')* p , l {r"A") (5.82) 

with wave functions as given in equations (5.79) and (5.80) and energy spectrum 
equations (5.78) and (5.81), respectively. 



5.3. Coulomb Potential and 1 /r -Potential in D Dimensions [15, 55, 91] 

In this Section we discuss the 1/r-potential in D dimensions. For D = 3 we have, of 
course, the Coulomb problem. The whole calculation is quite similar as in the previous 
cases, we just have to use a Kustaanheimo-Stiefel transformation in one dimension. 

The .D-dimensional path integral problem was first discussed by Chetouani and 
Hammann [15], whereas the radial problem by Ho and Inomata [55] and [91]. 

We start with the D-dimensional path integral with the singular potential V(r) = 
—QiQ2/r, where r = \x\ (x G R D ): 



x(t")=x" 



K(x",x';T) 



Vx(t) exp 



x(t')=x' 



m 

— x* + 



2 1 w dt 



(5.83) 



We denote the coupling by q\ x q 2 to emphasize the possibility that two "charges" q\ 
and q2 can interact which each other. As we have studied very explicitly in Section 
III. 3.1, the angular variables can be integrated out yielding 



K(x",x';T) = s[ D \{6"})S\ D) \{6'})Ki(r", r'; T), 



(D) 



(5.84) 



1=0 

where the radial kernel Ki(T) is given by 

r(t")=r" 

= (r'r") 1 ^ J £>r(t)exp 



r(t')=r' 
r(t")=r" 



III 



r 2 -h- 



2mr 2 



1 ■ q ^)dt 
r J 



+ 



dt 



r(t')=r' 



(5.85) 
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The ({0}) are the hyperspherical harmonics on the S D 1 -sphere. The correspond- 
ing Hamiltonian to the path integral (5.85) has the form 



H = 



h 2 ( d 2 d - 1 d 



^7 + 



+ h 



2m \ dr 2 ' r dr 
2m Pr + 



2 l(l + D-2) _ Q1Q2 

2mr 2 r 
QiQ2 



(5.86) 



2mr z 



with the momentum operators equation (3.45). We now perform the space-time 
transformation 

/,T s" = s{t"), r = F(u)=u 2 (5.87) 



8(t) 



4r(<r)' 



with / = 4k 2 , thus that F' 2 = f and we get the space-time formed Hamiltonian 



. 9 9 (2Z + D-2) 2 - \ 9 
H = r—Vu + ~ — -T ~ - 4Eu 2 - 4 qi q 2 



1 

2m 1 



with the momentum operator 



Pi. 



2mu 2 



2D -3 
i \ du 2u 



(5.88) 



(5.89) 



Here we have applied the results of space-time transformations for one-dimensional 
Hamiltonians. Thus we get the transformation formulas 



1 r°° 

K l (r",r';T) = -±- / dEe- lTE / h G l {r»y ] E) 
Znih J_ 00 



-oo 



/•OO 

G l (r",r';E) = 2i(u'u") 3 *- D / ds" e 4i ^ s "/ h Kt(u" , u'; s"), 

Jo 



(5.90) 



where the kernel Ki(s") is given by 
Ki[u ,u;s ) 

u{t")=u" 

= J Vu(s)fi 2 i+D-2[u 2 }exp 

u{t')=u' 



rn 



u 2 + AEu A ) ds 



n~n muJ 

V u u — — exp 

lnsmojT 



^(u' 2 + u> 2 ) cot uT 
Zi i n 



I-21+D-2 



muju'u" 



i h sin uT 



(5.91) 



with uj = ^—8E/m. The Green function we obtain in the usual way with the result 

Gi(r",r';E) 



(r'r") 



x exp 



2-d 2muj 



h 



F 

Jo 



ds" 



smws' 



4iQig 2 s" moj , „ „ 

— (r + r ) cot us 

n 2ni 



l 2l+D-2 



y/ r'r" 



ih sinus" 



I W? i D—l 9i 12 I rn~\ 

/ rn + — rT\/-2E ) 

(2Z + D-2)! 

8mE 



X W qi q 2 I m 7 I D- 
K V 2E''~ r 2 



> | M q 1 q 2 I m 1 i D- 
I K V 2E ' t_r 2 



8mE 



< 



(5.92) 
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Thus it is easy to display the wave functions and spectrum for the discrete contribution 
for the D-dimensional 1/r-problem [a = h 2 / (m\qiq2\) , n G N]: 



*N,i(r,{9}) = 



(N-l-1) 



2(N+^)(N + l + D-3)l\ \ a (N+2=Z) 



D/2 



2r 



x exp 



■ (21+D-2) 
'N-l-1 



2r 



a(N+^) 



a(N+^) 
SI D \{6}) 

(5.93a) 



(-1) 



N-l-l 



N + 



D-3 



1 



a(N + 1 + D -3)\(N - L-l)\ 



xr °^W N+£ -3 l+ n , 



2r 



D-3\ I I 



En 



mq\q\ 



For the continuous spectrum we get similarly 



(N = l,2,...) 



si D \{9}) 



(5.93b) 
(5.94) 



2n(2l + D - 2)1 



x exp 



(wH^ ( - 2i ^ <D,(W) (5.95) 



h 2 2 

with E p = -p^—- The complete Green function thus has the form 



J P 2m 

oo oo 

G(x",x';E) = X]E ex P iT 

AT=1 1=0 

°° poo 

+Y. ex p 

1=0 J ° 



mq\q\ 



2h 3 (N+^) 2 



nAr\{6'})* N Ar"Ae"}) 



inT £)^ (r, ' m) ^ (r// ' { ^ }) - (5 - 96) 



Finally we can state the following path integral identity [u = ^/—8E/m) : 



Vr(t)fj,\[r 2 ] exp 
2mu)\/r'r" 



i ivh 



/_ 



dE e" 



i TEh 



X 



00 ds" 



sincus' 



exp 



Mas" mu , n 



h 2\h 



[r 1 + r") cot us" 



muyr'r 
i h sin uis" 



2E T(2A + 1) 



76 



III.5.3 



Coulomb Potential and 1/r-Potential in D Dimensions 



x W q 



I rn \ 
"V 2E' 



8mE 



8mE 



ihj Q dE e~ iTE J E ^-L_^( r ')*„(r") + ^ 



(5.97) 



The discrete state wave functions and energy spectrum are given by (a = fi jraqxq^) 



*n(r) = 



1 



n! 



ar(n + 2A + 1) 
r 



2r 



a(n + A + |) 



x exp 



(-1)' 



n + X + \ 



a(n + A + \) 
1 

an!r(n + 2A + 1) 



2r 



L (2A) 



n+A+i A 



2r 



5 ' A Va(n + A + i) 



— 



ma 



2h\n + \ + i)2 



(5.98a) 

(5.98b) 
(5.99) 



For the continuous states we have similarly for the wave functions and the energy 
spectrum 



* p (r) 




rr(A+i-JL 



exp 



7T 



2ap 



2tt T(2A + 1) 

j-2 2 

with E p = -2^- • In particular this gives the identities 
m r(| + A - /j,) 



M±_ A (-2ipr) 



(5.100) 



ift/c T(1 + 2A) 

m \ -v n! 



E 



W M)A (2iA;r > )M At , A (2iA;r < ) 
1 



ft*; ^ T{n + 2A + 1) n + A + ± - \i 

x exp [ - k(r' + r")] (4/cVr") A+ ^i 2A) (2A;r')4 2A) (2fcr") 
1 m Z" 00 dp e^P 
2n hk J p — i \x 

— ^M ip , A (-2 i kr")M_ ip>A (2 i fcr'), 



x 



r 2 (i + 2A) 

or alternatively (and explicitly, a = (qi^/^) \/— m/2E) 
m r(| + A — a 



(5.101) 



2£ r(i + 2A) _Wq ' a 



8m£ 

) W 



8mE 
^T r < 



h 

00 dp 



1 /*00 7 

= ^^W'') + ?l / o E^E^'^^- (5 - 102) 
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Let us note that it is not difficult to analyze the Green function by taking the n th 
residuum in equation (5.92). Using the appropriate expansion for the T-function again 
we get 



1 



n=0 ^ «+^+A-2f 



2r^ 



2E, 

2r. 



m 



1 



2E n T(2A + 1) 



x W n+x+h ,^ a{n + j M n+x+h ^ a{n + j + regular terms (5.103) 



which gives using the relations of the Whittaker-functions with the Laguerre-polyno- 
mials the wave functions (5.92) and the energy spectrum (5.94). In particular we have 
for D = 1 



En 



1 



" 2r e- r ' an LP( — Y 



an° / an 
mZ 2 e 4 



n 



an J 



2h 2 n 2 



n G N 



l(i) = ^r 1-- e^Mi i(-2ipr). 
v 7 ^ V a P/ ap ' 2 



(5.104) 

(5.105) 
(5.106) 



Because, the domains x < and a; > are separated the wave-functions (5.104,5.106) 

are doubly degenerated. 

We can study several special cases (see [15]): 

i) We first consider D = 2 and get 



G{x",x';E) 



(r'r")"i [m~ g £(|Z| - ia^VS) 



2tt 




2£ 



00 



(2|/|)! 



9192 /m l; 
"H V 2E 'I 'I 



8m£ 



M i 



9192 /m 



,1*1 



8mE 



2m f 



sinhw 

oo 



exp 



9i 92 
ft 



2m v ' —2mE . ... 

-it ; (r +r Jcothw 



E 



h 



X 



m 



2tt ^ I ftsinhw / 

7 — rvi x * 



00 du 



7r/^ sinh-u 



cosh I z . ; vrr cos — 



h sinh it 



x exp 



h 



2m 



E 



■u — 



-2mE 



ir' + r") cosh-u 



(5.107) 



This is equivalent with equation (5.29). 
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ii) We consider D arbitrary. The Green function reads in the integral representation, 
where we expand in addition the hyperspherical harmonics in terms of Gegenbauer- 
polynomials 

nt ft i tti\ 2m / / f°° du 

G(x ,x;E) = -(rr ) . jf — exp 

oo 



9i 92 

ft 



2m \/—2mE . ... 

u - (r +r j coth-u 



E 



h 



X 



(2|^)£(g) ^ (12) 



2 v / -2mEVr 7 r" 
ft sinh it 



. (5.108) 



The /-summation can be performed by means of equation (5.74) and taking into 
account that we have for the Gegenbauer-polynomials 



(5.109) 



Thus we get 
G(x",x';E) 

2m 



= ^(r'r") — (4n) 



1-D 

2 



du 



sinh-u 



exp 



Qi 92 

ft 



2m \J—2mE . ... 

-it : (r +r ) coth-u 



ft 



x £ ( _ 1)1(2i + _ 2) E^zi> 



2=0 



x 2 -Pi [l + D — 2, -/; 



D — 1 



-; cos 



^,(1.2) 



2Z+d-2 



2 v / -2m£v / rV" 
ft sinh it 



2m/ 47rft \ 2 / r — V ( 

— — . v r'r" cos 

ft \^-2mEJ V 2 



2)X^ 



X 



(in 



o (sinhw) 2 



— exp 



9i 92 

ft 



2m V— 2mE . ... 

-u - (r + r Jcothit 



E 



h 



f2^/-2mE^fr T r" V (1 ' 2) 

X I d-3 cos 

2 y nsmhu 2 

Alternatively this can be written as 

G(x",x';E) 



(5.110) 



4nh \ 2 ( r'r" + x' ■ x" 



2m 

ft \V-2mEj 
_ 9i92 



du 



'o (sinh it) 



£+1 
2 



xexp 



ft V E 
m / 2-zrft \ 



2m \J—2mE . ... 

u ; (r +r Jcoth-u 



ft 



, 2V-2mE 
~ V ft sinh it 



r'r" + x' • x" 



l-D 
2 



q\q2 I 2m v ' —2mE 



du 



(sinhw) i 



x exp 



ft 



-u 



ft 



x coth it 



ft sinh it 



, (5.111) 
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where x = r' + r" and y = \x" — x'\. For D = 2 we recover, of course, equations 
(5.29,5.107) with I_i(z) = (2/tvz)? cosh 2. 
D = 1 yields with ipV* 2 ) -> 2tt: 



G(x",x';E) 

ft 7o 



sinh-u 



x exp 



ft 



u (x +x Jcothit 



E 



h 



h 



2E { h V 2E 



x W qi g 2 rrmi 1 1 \ h 

H V 2E ' 2 V V 



8m£ , ( ( . 

—9 T> I -M gl q 2 / m~~ 1 

^ / V~2B'2 



Furthermore one can show, by using the relation: 



2^/-2mEVx'x" 
h sinh u 



8mE 



Ti 



(5.112) 



d 

zdz 



1 v-\-m 



that the D-dimensional Green function for the 1/r-potential is connected with the 
(D — 2) -dimensional Green function by the relation 



G D (x\y; E) 



1 d 

2ny dy 



G D - 2 (x\y;E). 



(5.113) 



Therefore 



0-1 
2 



G D (x\y;E) = (3^) G Mv\E), 15 = 1,3,5,..., 



D-2 

2 



G D (x\y; E) = (^^) G 2 (x|y; E), D = 2,4,6, 



(5.114) 
(5.115) 
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5.4. Axially Symmetric Coulomb-Like Potential [7-13, 46-48, 88, 89] 

It is possible to consider the even more complicated potential problem, namely [46] 

,2 



to 



-ci 



q 



(x 2 + y 2 + z 2 ) + - 
z \Jx z + y z + z z 



l 



2mx 2 +y 2 2m ( x 2 + y 2) ^ x 2 + y 2 + z 2 



m,., 9 - 9 9 , q 2 bfi 2 c% 2 cos 9 

— (r 2 + r 2 9 2 + r 2 sin 2 96 2 ) + - = 9— 

2 r 2mr 2 sin 2 9 2m r 2 sin 2 9 



(5.116a) 



(5.116b) 



9w>9 -9\ ""1^9 9^9 2(/ 2 bf? C% 2 f] 2 — f 2 

- + ^ ^ + V 2 ) + T ^V^ 2 + 7^-^ " 7) ~PT~2 c2 L2I 2V 

2 2 t, z + t] z 2mt, z t] z m £ 77 (£ + 7/ ) 

(5.116c) 

Here I have displayed the corresponding classical Lagrangian in cartesian, polar and 
parabolic coordinates. It belongs to the class of potentials mentioned long ago by 
Makarov, Smorodinsky, Valiev and Winternitz [71] which are separable in a specific 
coordinate system and furthermore exactly solvable. This potential was discussed by 
Carpio-Bernido, Bernido and Inomata [10] starting from cartesian coordinates and 
using the four-dimensional realization of the Kustaanheimo-Stiefel transformation. 
However, this potential is separable in polar as well as parabolic coordinates and we 
shall give the explicit calculation for all of them. 



1) Cartesian coordinates 

In our line of reasoning we follow reference [10], however with some simplifications be- 
cause we have already discussed the realization of the four-dimensional Kustaanheimo- 
Stiefel in the discussion of the hydrogen atom. We consider the path integral repre- 
sentation corresponding to equation (5.116b) 



3N N-l - 

K(x",x';T) = lim ( -^V) 2 TT dx^dy^dz^ 
N^oo \2-K\eTiJ - LJ - / 



x exp 



. N 

-Y 



i=i 



^(A 2 x^+A 2 y^+A 2 z^)+e- r 
2e rw 



cbhf 



tchf 



Ai) 



2m rC?) 2 - z(J') 2 2m rC?) (rO') 2 - z^) 2 ) 



(5.117) 



Similarly to equation (5.45) we insert a factor "one" 



1 = lim 



TO 



N^oo \27Tieh 



v N 
2 



n 



d^ exp 



1 TO 

27n 



A 2 ^0') 



(5.118) 



81 



Important Examples 



We repeat the steps from equations (5.49) to (5.57), i.e. we realise the four-dimensional 
Kustaanheimo-Stiefel transformation on midpoints and arrive together with the 
transformation formulae 



1 r°° 

K(x",x';T) = ^^ J dEe- iTE / n G(x",x';E) 

/•OO 

G(x", x'; E) = i / ds" e Ai( ? s " ' h K{u" , s"), 
Jo 



(5.119) 



with the transformed path integral 



K(u",u'; s ) 



1 

4 r " 



^"K^u'^u^u^u^s") x K 2 (ulu' 3 ,ulu' 4 ;s"), (5.120) 



and the kernels Ki(s"), K 2 (s"), respectively, are given by 



K^u'l^ulu'^s") 

/.ff\ II (4-U\ It 

Ul(t )=U 1 U2(t )=u 2 



Ul(t')=M'i 

x exp 

Pl(t")=p'l 



Pui(s) y % 2 (s) 

U2(t') =,u 2 



^ Jo 



rn 



(u\ + "1) + 4£; K + "2) - 

Mt")=*i' 



2 , 2^ (6 + c)^ 



2 v_1 ' ' — v_1 ' 2m(u 2 1 +u 2 2 ) 



ds 



= f PiDpM J P0i(s) 



Pi(t')=pi 
x exp 



i 1 (*')=0i 



m ' -.2 i „2j^ i /iiti„2 t^2^"'" C 4 



1 -i 



2 -W + ^) + 4B^-V 2mp? 



(is 



1 oo 

* V e 1 



vi(<t>"-<t>'i) 



<ll\ // 



f,ll\ II 

Pl(t )=P1 



X 



Vpi(s) exp 



(5.121) 



-pl + AEpl-h 2 - 



Pl(t')=P' 1 



2mp\ 



ds 



2n i h sincus 



M<t>"-4>'i) 



v 1 = — oo 



x exp 



2 i n 



muipip'l 
i ft sinws" 



(5.122) 



Here we have introduced two-dimensional polar coordinates iti = picos^i, u 2 = 
pisin^i, separated the 0i and the pi path integration in the usual way, applied 
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the solution of the radial harmonic oscillator and have used the abbreviations u> = 
y/-8E/m, Ai = +V^ 2 + b + c. Similarly 



moo 



2tv i h sinus' 



E e' 



V2 = — OO 



V 2 (<f>2 -<t>2) 



x exp 



2^yP2 + P2 ) COt US 



i\ 2 7 



mup' 2 p% 



ih sinus 



77) (5-123) 



with us = P2 cos 02) U4 = P2 sin 02, A2 = +V v 2 + b — c. We identify variables [Cayley- 
Klein parameters, c.f.equation (5.59)] 



6 9 
Pi = Vr cos - p 2 = Vr sin - 



0i 



a + 



02 



a — 



(5.124) 



2 2 

(0 < < 7T, < < 2tt, < a < 4tt, = r"da"). Collecting factors we thus obtain 

2 



G(x'\x';E) = - / r/.s" . , „ 
7r J V 2 in, sin ws" 



exp 



4 i q 2 s" mu , „ 



— (r' + r") cot us" 



h 2\h 



x 



^ Ui^sinws" 2 2 / 2 V i 



— cos — cos — 

Usmus" 2 2 



(5.125) 

(Compare with equation (5.69)). Using the "addition theorem " [31, Vol.11, p. 99]: 

Z 

-J^,(z cos a cos (3)J u (z sin a sin (3) 
= (sinasin/3)^(cosQ;cos/3) M 

OI lN r(Ai + z/ + Z + l)r(z/ + Z + l) T 

x + , + 21 + 1) ^ r2( , +1)r( J ,;, +1) ^ m+iw 

x 2 i r i(-i, + + / + sin 2 a) 2 Fi(-Z, + + / + sin 2 /?). 

(5.126) 



0' . 0' 
sin — sin — 
2 2 



we get in the usual way by performing he s" integration 

G(x",x';E) 

v= — oo ^ 

vvj, ; x , 9 , , 1 , r(A 1 + A 2 + / + i)r(A 1 + / + i) 
t^ (Al + 2 + + !) /-r^ + i^ + z + i) 



0' 0" 
cos — cos — 
2 2 



z=o 
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x 2 Fi ( - I, Ai + A 2 + I + 1; Ai + 1; sin 2 °— J 2 fJ - /, A x + A 2 + / + 1; Ai + 1; sin 2 °— 



x 



00 cte" 



\/r'r"U Jo sin cus" 
x exp 



— -(r +r") cot ws" 



mu)\Jr'r" 

-'Ai+A 2 +2Z+l 7^-7 



i ft, sin cus" 

2tt ^ e f-' "2 r(Ai + / + i)r(A 2 + / + 





f . 0' 


6"\ Xl i 


' 0' 






sin yJ ( 


cos — COS 




V Sm 2 




2 



>(Ai,A 2 ) / n'\ 73(^1,^2) , 



X 



m r[J(Ai + A 2 ) + / + l-p] 



r'r" V r(Ai + A 2 + 2/ + 2) 

x Wp |tH i±xi±Aa (-2iA;r>)M p|tH i+a 1+ a 2 (-2 i fcr<). 



(5.127) 



with the abbreviations p = q 2 y — m/2Eh 2 , k = y/ +2mE/h. This is the result of 
reference [10]. The energy-levels are determined by the poles of the T-function and 
are given by 

4 

mq 



E n — 



2ft 2 (n + Z + l + Ai±A2)2' 



n e N . 



(5.128) 



The corresponding wave functions will be calculated below, where we separate the 
radial from the angular path integrations directly. 

2) Polar coordinates 

We consider the path integral corresponding the Lagrangian (5.116a): 



K(x",x';T) 

r(t")=r" 



e{t")=e' 



4>{t")=4>" 



= J Vr(t) J sm0V0(t) J V(j){t) 



r(t')=r' 0(t')=0' 
-t" 



x exp 



i 

ft .It' 



<t>(t')=4>> 
?n (r 2 + r 2 2 + r 2 sin 2 0</> 2 ) 



ch 2 cos 9 U 2 



dt 



r 2mr 2 sin 2 9 2m r 2 sin 2 9 8m \ sin 2 9 . 

1 oo 

= — Y eW'-*') K v (r",r',6",e';T) (5.129) 
2n L — ' 

v= — oo 

with the kernel K„(T), which in turn is also separated: 
tf„(r",r',0",0';T) 
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r{t")=r' 



e{t")=e' 



= (r'r" sin 0' sin 0") - ^ J Vr{t) J V9{t) 



r(t')=r' 



9{t')=8' 



x exp 



m 
i 



(r 2 +r 2 9 2 ) - Ti- 



,v 2 + b + ccos9 - \ q 2 h 



2mr 2 sin 2 9 



+ — + 



r 8mr 2 



dt 



(r'r" sin 6' sinfl")" 5 

r(t")=r" 6{t")=6" 



X 



rVr(t) 



V6(t) exp { - 

Ti 1+/ 



r(t')=r' 



0(t')=6' 



m 9 Am o / 9 
-r + — r 2 - 



~z/ 2 + & + c-± „z/ 2 + b- c-\ q 2 h 
- h — » : o - h — s ttzt- + — + 



8mr 2 sin 2 9 



8mr 2 cos 2 9 



8mr 2 



dt 



Y,^\ XlM \0')^\ XlM) {0")K lv {r",r';T), 



(5.130) 



1=0 



the (Poschl- Teller) wave functions 



^(A 2 ,Ai) 



(0) = 



"A 1 + A 2 + 2/ + l /!r(A! + A 2 + / + l) 
2 i+a 1+ a 2 r(X 1 + l + l)T(X 2 + l + l) 

x (1 - cos 0)^(1 + cos6>)^p/ Ai ' A2) (cos6>), 



(5.131) 



(Ai )2 as before) and the radial path integral 



KUr",r>;T) 




m .9 q 2 

—r 2 + - h 2 

2 r 



2mr 2 



dt 



r[t )=r 

^7, J Vr(t)exp 

r(t>)=r> 

(5.132) 

[A = 2(^1 + A2 + 2/ + 1]] . This path integral is nothing but a Coulomb potential path 
integral, however with generalized angular momentum. Thus we just analytically 
continue the known result in / and obtain 



G> (r", r'; E) = -L ^ |^ V|(-2 i krJM,^ (-2 i *r<). 

(5.133) 

The energy spectrum was already stated in equation (5.128) and the corresponding 
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wave functions are 

* A (r,M) = ^*! Al ' A2) (0) 



2tt 
2 



(n + A + i) 2 



n! 



a 3 (n + A + |)r(n + 2A + 2) 



x exp 




rr(A + i + -M 



x < 



2 ap 

2tt rr(2A + 2) 
and a denotes the Bohr radius. 

3) Parabolic coordinates 



2r 



a(n + A+ |) 



W2A+1) 



a(n + A + |) 
2r 

o(n + A + |) 



(5.134) 



6XP (i) M ^^ ( " 2iPr) 



The potential described in equation (5.116c) is also separable in parabolic coordinates, 
in the operator approach as well in the path integral formalism. We consider the path 
integral formulation corresponding to the Lagrangian (5.116c): 



K(x , x ; T) = K{£ , £ , 77 , 77 , , ; T) 

€(*")=€" v(t")=v" W')=<t>" 



Vcj>(t) 



4>{t>)=4>' 



1 



xexp^- 



y(e 2 + ^)(^ + r) 2 ) + ^V^ 2 



2 9 2 eft 2 r/ 2 -^ 



dt 



£ 2 + r/ 2 2m£ t 2 rf m £ 2 r/ 2 (£ 2 + rf) 

U= — QO 

with the kernel K V (T) 

€(*")=*" •»(*")=»/" 

^(e // ,e / ,r7 // ,r 7 ';T) = (erw / )- i / / ^ 2 +?? 2 ) 



*»(*')=»»' 



xexp^- 



+ 



2q s 



ft 



o6 + Z/ 2 -i c ^2 r] 2_£2 



£ 2 + r/ 2 2m£ t 2 rf m £ 2 r/ 2 (£ 2 + ?7 2 ) 

This path integral is now tracked by a time-transformation according to 

da 



)+V 2 (ay 



(5.135) 



dt } (5.136) 



(5.137) 
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therefore 



1 r°° 

K V {Z'\ 17", 1/; T ) = ^-J dE e- iET ' h G„(£", 17", 1/; £?) 

/■OO 

G„(£",£',r/',r/;i<;) = i / ^"e 21 ^^^,^^^,") 

^0 



(5.138) 



with (note that this is a two-dimensional time-transformation and the prefactor is 
"one"): 

= ^i(£",£V) x k 2 {rf\rf;a"). (5.139) 

This gives for the kernels K lj2 (s // ) 



*(*")=€" 



^ Al [e 2 ]^e(s)ex P 



(is 



iftsincus 



77 ex P 



mw , .,2 .,2 N // 

2ifc^ + 



i ft sinks' 



(5.140) 



with u = y/-2E/m, Ai = + V^ 2 + 6 + c. Similarly for K 2 {s") 

r)(t")=r)" 



H\ 2 [rf] exp 



v(t')=v' 



nl 



rn 



ff + Erf 



ds 



i ft sinews 



77 exp 



— — (77 +77 cotws 
2 1 ft 



mojrj'rj" 
i ft sinws" 



(A2 = +V^ 2 + 6 — c). Integrating over s" yields 

1 ' ' ' 2-^2.^ E N -E 



V = — O0 m ,Tl2=0 

with the energy spectrum (note to = q 2 /UN) 



(5.141) 



(5.142) 



En = 



mq 



2ft 2 iV 2 ' 

The wave functions have the form 



AT = m + n 2 + 1 + 



Ai + A 2 



(5.143) 



Ai 



2tt 

A 2 



2 2m!n 2 ! 
a 2 iV 3 ' r(m +Ai + i)r(n 2 + A 2 + l) 



£ 2 + ?7 2 
2aiV 
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For the continuous spectrum we obtain 



G <-vv;£)=a t r * r d( ^M^i^Mi (5 , 45) 



with 



B P = ^ (5.146) 



and the wave functions 

e i! ^ 1 L— 5- r r ~ >\< L 2~~ r 2 



(27r)i ^r(l + A!)r(l + A 2 ) 

_,_ 2 ) _ - 

2 ^ s 1 ap 2 2 » ap -" 2 



x M_ i f<r+ _M m (- i^ 2 )^_ i (c -^) h2 (- i PV 2 )- (5-147) 



Finally we state the relation of the Green function in parabolic coordinates. We have 

2 oo 



(\ z oo 
' !^ = — OO 

/•°° cte" _ / mw^" \ _ f wwrf-q" 

X ' * J >1 ^~T^ « /A 2 



x exp 



o sin 2 us" 1 V i ^ sm UJS " J 2 \ihsmus" 
2iq 2 „ mu /2 t ,i ,2 , ,2 



, s "-^(r+r+v z +^)coto; S " 



(5.148) 



Note the similarity to the u±, U2, u^, approach. Actually the Kustaanheimo-Stiefel 
approach in cartesian coordinates produces a separation in parabolic coordinates. We 
switch bach to polar coordinates by means of 

£ 2 = r + z = r(l + cos 9) = 2r cos 2 - 

2 (5.149) 
rf = r — z = r(l — cos 6) = 2r sin 2 - 

Using again the addition theorem (5.126) with a rescaling s" — > 2s", cj — > w/2 so that 
to = y/—8E/m we recover the Green function (5.125). 

Let us finally note that even more complicated Coulomb-like potentials can be 
exactly solved [47, 48] which are, however, only separable in parabolic coordinates. 
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